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ABSTRACT 



In the measurement of microwave circuits and devices, the data reduction of two- 
port measurement data obtained at the external terminals of microstrip-coupled devices, 
with correction for measurement error, can be obtained by de-embedment procedures 
with use of the HP 8409B vector network analyzer and a Fortran program. 
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I. INTRODUCTION 



A. BACKGROUND 

Microwave circuit elements, such as semiconductor devices, cannot be connected 
directly to the terminals of a microwave network, analyzer for measurement of scattering 
or impedance coefficients but must be coupled through a chain of hardware elements 
comprising connector-adapters, microstrip sections, and bondwire connections. The 
basic problem under consideration is the removal of the effects of intervening hardware 
from measurements made at the external terminals of an embedded device or circuit 
component to obtain the intrinsic device characteristics. Accurate measurement data 
on a device is essential for successful circuit design. 

Earlier work on device de-embedment has been coupled with undesirable restrictions 
such as the assumption of physical symmetry of hardware, inadequate provision for re- 
alistic measurement standards, or neglect of statistical treatment of data for improve- 
ment of accuracy. Consequently, a need exists for a practical treatment of this 
measurement problem which can be implemented in typical laboratory circumstances, 
with acceptable accuracy of results. 

In order to measure the scattering parameters of a device, the device is often 
mounted in a fixture. This future has a set of scattering parameters associated with it 
so that the parameters measured at the terminals of the future are the scattering pa- 
rameters of the device transformed by the embedding network of the fixture. 

Thus, to obtain the scattering parameters of the device itself, the device must be 
"de-embedded" from the future. Techniques for using a network analyzer (N.A.) to de- 
embed the device from the test fixture have been discussed earlier (Refs. 1,2,10,11]. 

As long as the device interface connection is coaxial line or waveguide, de- 
embedding is relatively straightforu'ard because standards are available. But if the device 
is constructed without a coaxial or waveguide interface, as in the case of microstrip or 
coplanar line systems, de-embedding becomes more difficult. 

Extracting the device parameters from the overall measured parameters is strongly 
dependent on the ability to characterize the futuring. Figure 1 is the block diagram and 
signal flow graph of a typical N.A.- futuring setup. The problem is to find the best 
characterization networks A and B so as to remove their contributions, yielding the de- 
vice parameters, S. 
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Figure 1. Block Diagram and Signal Flow Graph of a Typical Fixturing Set-up 



B. SCOPE 



The scope of the proposed research will be to develop a working procedure for ear- 
thing out the removal of embedded device scattering matrices from network analyzer 
data. The procedure developed will be supported by experimental measurements on 
typical microwave circuits fabricated on planar substrates. 
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II. CHARACTERIZATION OF THE FIXTURES 



A. CHARACTERIZATION CHOICES 

Several papers present characterization techniques that are in use at various labo- 
ratories. There are four main choices of N.A.-fixturing characterization as shown Fig. 
2 [Ref. 1]. 




Figure 2. The Four Main Choices of N.A.-Fixturing Characterization. 



Method (1) combines a variety of measurements, including the transmission between 
fLXtures connected back-to-back and reflection measurements on the fixtures terminated 
in microstrip stubs. In this way, data is assembled to enable the scattering coefficients 
of the embedded device to be converted to the device scattering matrix. Method (1) 
therefore contains features of both methods (3) and (4) below, and represents also an 
early e.xample of the use of microstrip stubs as known reference standards. 

In method (2), the measured S-parameters are fit to an assumed circuit model to 
calculate the device parameters. 
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Method (3) has the advantage that only three simple test pieces in the device me- 
dium are required. More details are given in Chapter V, part B. 

From the above choices, method (4) was selected for investigation. In this method, 
three or more known reflection coefficients are measured at each port. The reference 
load reflection coefficients are derived from independently measured or calculated non- 
zero values. Reference 2 used a chip capacitor, wire bonded to microstrip on alumina, 
with capacitance measured at low frequency (LF). The overall reflection coefficient was 
then measured versus frequency. To obtain multiple calibration measurements, a piece 
was scribed off the chip, and the measurement cycle repeated twice more at each port. 

In other cases, the reference load was a high-Q varactor chip, again measured as a 
function of three values of reverse bias at LF. The chip was then wire-bonded to the 
microstrip conductor and the overall reflection coefficients of the fixture-chip combina- 
tion versus frequency were measured at each value of reverse bias. 

In the present work, the reference load is a section of the microstrip line itself. Three 
different lengths of microstrip line yield the three different reflection coefficients for each 
port by separate measurements. 

B. THE TEST FIXTURE 

The S-parameters of the test fixture can be determined by measuring a set of known 
impedances in place of the embedded device. This process is known as "unterminating." 
Once the fuxture's S-parameters have been determined, the device is measured in the 
calibrated test fixture, and the device parameters are "de-embedded" from the data to 
give an accurate picture of device performance. 

The signal flow diagram and block diagram in Fig. 3 defme the input and output 
characteristics of the test fixture. 

Three sets of S-parameters must be determined when an unknown device is charac- 
terized as a two-port network mounted in a test fixture. First, the S-parameters of the 
input (A) and the output (B) sides of the fixture must be evaluated before the S- 
parameters of the device itself can be deduced. If the test fixtures are known to be sym- 
metrical, then only one measurement is needed to determine the S-parameters of the two 
sides of the fixture. 

In the present case, however, unsymmetrical fixtures are chosen, requiring six dif- 
ferent measurements and calculations for both fixtures independently. To solve this 
problem, the input half of the fixture is first terminated at the device plane with an ar- 
bitrary load of reflection coefficient, F^^, followed by other loads with different reflection 
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Figure 3. The Flow Diagram and Block Diagram of the Fixture, 
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coefTicienls. Then the output half of the fixture must also be considered with same 
method. 



a1 



a2 







Figure 4. One Port De-Embedment 



For the one port case, as shown above, the reflection coefficient is the data needed 
since the impedance can be calculated from it; 

ai) 

From the definition of the scattering matrix, the input reflection coeflicient is: 




«i 

where, for load r^, 

= F [b2 (-•^) 

and 



— '^ 21^1 ^ ^ 22^2 



(2.4) 



— S 2 ]£ Ii -b 522(^1^2) 



1 



Thus, 



b, = S- 



21 



1 - s„r 



22 * L 



Using the above, we get: 

■^12‘^21^ L 



^1 1 ~ 1 T 






(2.4a) 



(2.5) 



Here, S,, can be found by measuring a matched load in the fixture (F^, = 0). If we 
assume that this measurement gives F„ = then we have S„ = F^^. Using the result 
in Eq. (2.5), 

(H/i ~ ^.V/l)(1 ~ ■^22^z3 = •^12-^21^'l (2-6) 



From this. 



■^i2'^2i^~i: + ■^22^ jn/i ~ — (Hn “ I'wl) 



(2.7) 



Two more measurements with known loads will provide enough data to determine 
and (SjjSj,). Suppose we measure with an open-circuited load (Fi. = +1). and a 
short-circuited load (F^ = —1). Then, 

■^1 2-^21 + •^22(^~OC ~ ~ ^ML 

~ ^22i^ SC~^ ml) = fsc “ r'wL 
These can be solved for 

f'oc + ^" 5 ^ — 2 F^£^ 



S 22 — 



•^12-^21 — 2 



f'oc “ ^5C 

(Toe ~ J ML ~ ^5c) 

foe “ f'se 



can be written from Eq. (2.7) 

(n„-F,v,j 



F, = 



^ S,2S2i + S 22 (F i„ — F 



(for 


rz, = +i) 


(2.8) 


(for 


rz, = -i) 


(2.9) 






(2.10) 






(2.11) 


1 for de-embedding t 


i general unknown 






(2.12) 
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But in the actual case, it is very dilTicult to construct the matched load and ideal 
shorted and open terminations required by this procedure. Therefore, we rearrange Eq. 
(2.12) to matri.x form for three different known loads and their reflection coefficients. 
From Eq. (2.5), 




•^11 -^1 2-^21 

1 - s.2rL 



_ S, 1 •^1 1-^22 + ■^12‘^21 

J_L 

1 - S22fL 



Rearranging this gives; 



(2.13) 



7^ = -r^ + - iSuS22 - SnS2i) 

1 £. 1 L 

(2.14) 

^11 

= -p^ + r,,S22-A 

i L 



where, 

^ = -^1 1-^22 ~ ■^12'^21- (2-15) 

This equation then involves the three unknowns S,„ S^, and A . These three unknowns 
can be determined by making three independent measurements off,, with three different 
known loads. In matrix form, the problem can be expressed as: 



1 


i/ii 


-1 




•^11 




ln\ 


Tu 






Tl. 


1 




-1 




S22 




f~(Vi2 


r^2 


* in 2 






Tz.2 


1 












nn3 




H/is 


-1 




A 





To find S„, S 22 and A, solve: 



(2.16) 



9 



-1 



( 2 . 17 ) 



■^11 

$22 

A 




T/nl 

r m2 



-1 



1" in\ 

ru 

1" ln2 

ri2 

r m2 



f ill 1~ L2 and r^3 are three difierent calculated reflection coefficients at the terminals 
of the three different loads. Also, r„,, and r,„3 are three reflection coefficients 
measured from the N.A. with respect to the three different loads. Now the problem re- 
quires accurate values of Fi,, F^j, F^j and F,„i, F,„2, F„j. We consider in the next chapter 
the computation of F^„ F^^ and F^j for reference microstrip line loads. 
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III. MICROSTRIP LINE FOR THE FIXTURE 



A. MICROSTRIP LINE 

The definiiion of an effective relative dielectric constant, e,,, allows us to treat 
microstrip line like a TEM system. Relations between line and freespace parameters 
then hold (which are typical of TE.M systems). With replacing the conventional rel- 
ative dielectric constant, 



V ^re 



'’ms 

c 



^Oair 



V ^re 

I 



re 



V C 



^ms ^ms 



(3.1a) 

(3.16) 

(3.1c) 

(3.1^ 



where and v„, are the wavelength and velocity of the waves on the microstrip line and 
/to and c are the free space values at the frequency of operation. 

In Eqs. (3.1), L„, and C„, are the inductance and capacitance per unit length of the 
microstrip line. Zo is its characteristic impedance, and Zo,„ is the value Zq would have if 
the dielectric layer were replaced by air. To the lowest order of approximation, the value 
of £„ would be the average of the relative dielectric constant of the substrate and that 
of air, namely: 



£ 



re 



1 + £,. 
2 



(3.2) 



More exact values have been obtained, largely through solutions of the quasistatic 
model, corrected to match measured data. 
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Figure 5. Microstrip Line Dimensions. 



In case of lines with line width ratio « s -^ = 0.1 to 10, 

h 



£re(0) = 



e,. + 1 Cr ~ 1 



Ly 



1 + 



12.35 



+ 0.06(1 



(3.3) 



If microstrip were a true TEM system, the velocity of waves on the line, the 
impedance and other line parameters would be independent of frequency. The 
quasistatic or zero frequency values of the parameters are good approximations at fre- 
quencies up to the order of a few gigaherz. When the wavelength becomes comparable 
with the cross sectional dimensions of the line, and the line parameters are no longer 
frequency independent. Therefore, dispersion of effective relative dielectric constant is 



gr - gre(O) 

2^^ 2 
^ Zo ^ 



(3.3a) 



where h is the substrate height, f is frequency and is the permeability of free space. 

B. MICROSTRIP STUB-RESONANCE EVALUATION OF DIELECTRIC 
CONSTANT 

It was necessary to determine the substrate dielectric constant e, of the microstrip 
line substrates used for measurements in the present work. 
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The steps are: 

1. Vleasure the resonant frequency of a quarter wave length stub. (Substrate height 
is h= 1.5 mm, line width ratio u = 3.0 (where u = w / h)), 



4 Different microstrip stubs, with: 

— Length = One quarter wave length 

— Substrate height = 1.5 mm 

— w / h (Line width ratio) = 3.0 




Figure 6. Measurement Setup for Dielectric Constant Measurement 



2. Measure the value of Sjj. transmission coefficient, with respect to frequency using 
the Ni.A., 

3. Find the frequency {Q at which ma.ximum dB transmission loss due to resonance 
occurs, as shown in Fig. 7, 




Figure 7. Determination of the Max. dB Loss Point. 
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4. Calculate £„ from Eq. (3.1a). It is normally required to account for the fringing 
capacitance at an open microstrip termination. Computation by standard methods 
[Ref 3: p. 138] showed that the open-end length correction was less than 1.0% for 
stub lengths greater than 50 mm. This had negligible effect on the calculated £„ 
and was neglected. 

^•0 ~ ^re (^•‘^) 



™ /£ 

V ‘re 

therefore, 



) = — ^ 
/Ov ^ 



OV ^re 



where 



(3.5) 



(3.6) 



X 



ms 



= 4/. 



From the equations above, 



(3.7) 

Through the several measurements as shown in Table 1, the average e„ ob- 
tained by calculation has the value of 2.833., 



Table 1. THE MEASUREMENT OF EFFECTIVE RELATIVE 
DIELECTRIC CONSTANT 



Length (mm) 


UGHz) 




98.0 


0.452 


2.867 


80.0 


0.564 


2.763 


60.0 


0.746 


2.808 


40.0 


1.102 


2.895 



5. From stub resonance data, the corresponding is calculated. Eq. (3.3) was in- 
verted by trial and error to find e, = 3.516 from the experimental t„. 

C. CHARACTERISTIC IMPEDANCE OF MICROSTRIP LINE 

Expressions for the characteristic impedance of microstrip lines have been developed 
through numerous successive approximations, based principally on TEM models having 
the specific microstrip line cross section, with correction terms added to obtain agree- 
ment with experimental observation. 
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(3.8) 




(0.99^+ 1.795)(0.994f^ + 5.80) - 
u(u+ 1.268) 



where the synthesis formula for linevvidth ratio u = w;h is: 




^ ^/.J(0.647 + 



0.4116 



£ 



)+ 1.321(1 +-f) 

‘■r 



(3.9) 



with 




(3.10) 



From Eqs. 3.3, 3.6 and 3.8, one can easily get the characteristic impedance for each 
microstrip section. For example, for the given characteristic impedance of 50 Q , the 
microstrip line width ratio is u = 2.283 and the width is 3.43 mm. In the construction 
of microstrip stub loads, it is necessary' to choose high characteristic impedance for ac- 
curate measurement, since high impedance leads to narrow' linew'idth and a well-defined 
point of attachment of the stub. Construction problems arise in making microstrip line 
with very small width. For example, if we want to make a microstrip line stub of 90 Q 
characteristic impedance, the required width is about 1.2 mm for given h = 1.5 mm. It 
is difficult to handle the adhesive foil material used for the microstrip itself For this 
reason, a compromise value was chosen, using 2.0 mm width of microstrip stub for given 
substrate height h = 1.5 mm. Then, 

h = l.5mm 

w = 2.0mm 

also £/ = -r-= 1-33. 

h 

Therefore, 

Zo^ = 68. 2Q, 

which was chosen as a w'orking value for stub impedance. 
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MICROSTRIP UNE 



CHARACTERISTIC IMPEDANCE : 50 Ohm 
WIDTH : 3.43 mm 

OPEN STUB 

CHARACTERISTIC IMPEDANCE : 68.2 Ohm 
WIDTH : 2.00 mm 




Figure 8. Experimental Setup for Characterization of the Fixture 



D. STUB LOAD IMPEDANCE AND REFLECTION COEFFICIENT 

Consider the transmission-line impedance- transformation law; 

_ r cos /?/ sin /?/ ] 

|_ Zqj cos /?/ -f sin /J/ J ^ ^ ^ 

From Eq. (3.11), when Z^, goes to oo (open ended), the input impedance of an open • 
ended stub becomes: 
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(3.12) 



^In 




COS /?/ 

j sin pi 



= -JZqs cotiPl) 



using 




V 

f 

_ C 

f\! ^re 



Therefore, 




(3.13) 



(3.14) 



(3.15) 



Once the Z^, is chosen, we can determine the load impedance of three different loads. 
For example, given Zo, = 68.2 ohms and microstrip line lengths of 4 cm, 3 cm, 2 cm, the 
load impedances are tabulated as shown Table 2. 



Table 2. THE LOAD IMPEDANCE OF REFERENCE STUB LOADS 



Frequency(GHz) 


L = 4 cm (Zii) 


L = 3 cm (Z^i) 


L = 2 cm (Zij) 


1.00 


-i 11.06 


-j 38.44 


-j 80.15 


2.00 


j 204.8 


j 41.28 


-j 11.06 


3.00 


-j 35.69 


-j 1230 


j 41.28 


4.00 


j 91.05 


-j 35.69 


j 204.81 


5.00 


-i 70.73 


j 44.23 


-j 168.99 
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The three different reflection coefficients of the given arbitary loads are; 



r^i = 



~ •^0 



(3.16) 



where i = 1, 2. 3. Therefore, from the given stub impedances and Zo = 50Q, the re- 
flection coefficients are tabulated in Table 3. 



Table 3. THE STUB REFLECTION COEFFICIENTS 



Freq.(GHz) 


L = 4 cm (Ft,) 


L = 3 cm (Fj^,) 


L = 2 cm (Fj^j) 


1.00 


-0.91 - jO.42 


-0.26 - jO.97 


0.44 - j0.90 


2.00 


0.89 + j0.46 


-0.19 + j0.98 


-0.91 -jO.42 


3.00 


-0.32 - jO.95 


1.00 - j0.05 


-0.19 + j0.98 


4.00 


0.54-(-j0.84 


-0.32 - jO.95 


0.89 + j0.46 


5.00 


0.33 - jO.94 


-0.12 + J0.99 


0.84 - jO.54 
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IV. APPLICATION OF THE DE-EMBEDDING PROCEDURE 



A. MEASUREMENT OF THE EMBEDDING FIXTURE 

The reflection coefficients required are measured by the HP 8409B vector network- 
analyzer. One can find the reflection coefficient, f„„ of the stub terminated fixtures by 
measurement of the scattering matrix element Si, using the network-analyzer. The out- 
puts of the N.A. represent the magnitude (dB loss) and phase angle (degrees) of S,,. 
Therefore, the outputs of the N.A. must be changed to complex notation using the 
simple modification: 

~^ndB -^UdB 

n/i = (10) 20 cos( arg Sii) -l-y(lO) 20 sin(argS„). (4.1) 

The outputs of the N.A. and f*, are obtained for each load on one fixture (A) as shown 
in Figs. 9-a to 9-c. When we consider the other fixture (B), all steps are exactly same 
as for fixture A except the fixture shape itself, shown in Figs. 10-a to 10-c. Tables for 
fixture A and fixture B data are available in Appendix A. 

B. S- PARAMETERS OF THE FIXTURE 

To consider the de-embedding calculation, we return to the original matrix equation. 
The success of this characterization technique depends heavily on the measurement of 
the load impedances used in the unterminating process for the fixtures. Once an input 
reflection coefficient has been measured with three different impedance standards, the 
S-parameters of the test fixture can be computed using the matrix introduced earlier. 
A word of caution must be offered about the convergence of this apparently simple 
problem. 

The matrix that must be inverted consists of elements very near unity, a situation 
often used to illustrate some of the numerical inversion methods. In this case, the 
Gauss-Jordan method with full pivoting is recommended. Pivoting improves the stabil- 
ity of the numerical solution by interchanging rows and columns in the Gaussian elim- 
ination. In the notation of Eq. (2.17), the scattering parameter S,j is changed to A^j and 
for fixture A and B (respectively): 
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Figure 9. Measurement of the Load for Fixture A 
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Figure 10. Measurement of the Load for Fixture B 
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^11 
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HaiI 


-1 




Tli 






1 




-1 


^22 


— 


r^2 


1 m2 






1 






A 






f/AiS 


-1 



where 

A = A]iA22 — ^12"^21- 



r., 

^ ini 

r^2 

r/«3 



(4.2) 



(4.3) 



The computed values for .4,, and A 22 , combined with a knowledge of A , can be used to 
determine the A 21 product. Note that the fixture is reciprocal, so A^ = An and for the 
other fixture, = B^. Thus, the magnitudes of these terms can be calculated by taking 
the square root of the product term in the expression for A . 

The question of phase is left open, but is not essential. Since the A^i, 5,2 and 
Z? 2 i terms always appear in product form, it does not matter which of the two possible 
phase angles one chooses, as long as the phase angles for /4j, and B^ are approximately 
the same. All of these calculations are shown in Appendix E using the Fortran program 
showm. This program combines the measurement and calculated values, that is, the 
outputs of the N.A. are used in the calculation as a form of datafile. The resulting S- 
parameters of these fixtures are shown in Appendix B. 

C. DE-EMBEDDING PROCEDURE 

1. Background 

When a two-port network is connected between fixture A and B, the signal flow 
graph for the system may be drawn as shown in Fig. 11. Parameters AB and BA are 
mainly contributed by the finite isolation of the internal switches of the network ana- 
lyzer. The terms 5n, Sji, Su, and S 22 are the scattering parameters of the device under 
test. The fixture parameters. An, /lj„ and A 22 , can be determined by reflection 
measurements at port 1, as discussed earlier. Parameters 5jj, B^ and can be eval- 
uated by similar measurements at port 2. Two port transmission measurements with the 
device under test replaced by a through-line are used for evaluating parameters AB and 
BA, and products /4j,52, and A^B^- When a through-line of length /, having zero re- 
flections at both ends and a transmission T = c- "', is connected between the two ports, 
the signal-flow graph of Fig. 11 can be simplified to the one shown in Fig. 12. 
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Figure 11. Signal Flow Graph for Measurements on a Two - Port Network 



The measured transmission T„ can be expressed as 



Tm = + 



1 — 



( 4 . 4 ) 



When the error parameters A22 and B„ have been determined by the reflection measure- 
ments, two transmission measurements must be performed for the determination of BA 
and A2,B2,. The value of BA can be found by making a transmission measurement w'hen 
no connection exists between the ports (that is T = 0). The parameter /l2i^2i can be 
determined from the measurement of a through - line using 



where T = j/? being the propagation constant of the through-line. The error pa- 

rameters AB and A^2B^2 can be determined by making a similar measurement. 

Now that the scattering characteristics of the embedding networks have been 
measured, the scattering parameters M^^, M22, T/,2 and iV/21 can be corrected for the em- 
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Figure 12. Signal Flow Graph with a Through - Line 

bedding networks to obtain Sn, S^ 2 , 5,2 and Sj,. Elements M,j represent the scattering 
parameters of both fixtures and device, namely: 



,V/n = /4,i + 



Mj, = BA + 



M 22 = B 22 + 



iV/,2 = ^5 + 






D 



■^ 21 -^ 2 l -^21 

D 



•^22^21^12(1 — ■Si 1'»422) + (•S^21‘^12^22^12^2l) 



D 



■^ 12 -^ 12^12 

D 



(4.6) 

(4.7) 

(4.8) 

(4.9) 



where 



D — \ — S,,/l 22 “ •^22^11 “ -^l 2 •^21 '^22^11 •^ 22 ‘^ll'^ 22^11 



(4.10) 
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These equations must be inverted to extract the scattering parameters S^. An explicit 
solution for the scattering parameters of the device under test has been presented in Ref 
(4| : p. 313 . The solution is as follows: 

G(1 + 

^11= Y (4.11) 



s 

•^12 - Y 


(4.12) 


•>21 - Y 


(4.13) 


H{\ + A^^G)-A^^EF 
^22 - Y 


(4.14) 


where 




iV = (1 + A^^G){\ + 5„//) - A^^B,,EF 


(4.15) 


{ytn-AB) 
(^12^1 2 ) 


(4.16) 


(-^^21 - BA) 
{AM 


(4.17) 


("^12^21) 


(4.18) 


Y (-^^22 - ^22) 

{BnB2^) 


(4.19) 



In the present measurement, inclusion of the leakage paths BA and AB is not necessary 
because of the one-port calibration procedure used. Therefore, Eqs. (4.16, 4.17) are: 



Mn 

•^ 12^12 



(4.16a) 



F = 



^2:^2: 



(4.17a) 



Some variations of the above calibration procedure are available in Ref [4J : p. 315 . 
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The (Je-embcdding procedure used in this work is based on Oqs. (4,11 to 4.15) 
and Eqs. (4.16a, 4.17a, 4.18 and 4.19) which use the 9 key parameters; 

•^1I> ■^I2> ‘^21> •^22> 

and the foregoing derivations. The complete block diagram is shown in Tig. 14. If we 
know A,j and D,j, the scattering parameters of fixtures A and B, and the parameters of 
both fixtures plus device, M,j, then the scattering parameters of the device itself are ob- 
tained by using Eqs. 4.11 to 4.14. In order to test the de-embedding procedure, consider 
the set up of Fig. 13. 



NETWORK 

ANALYZER 



FIXTURE + DEVICE 



FIXTURE A FIXTURE 8 




Figure 13. Device and Fixture Assembly for the 51 Ohm Resistor 



It is necessary to choose a device that has known scattering parameters; then it will be 
possible to compare the output of the de-embedding experiment with the theoretical 
parameters of device itself. The initial test object is a series-connected 51 Ohm resistor 
which has a small inductance so that the impedance of the 51 Ohm device is: 



Z = R + j(oL 
=5=51 +J2n/L 



(4.20) 



where we may assume the parasitic inductance associated with the resistor to be 
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Fixture B 




Figure 14. The Block Diagram of the De- Embedding Procedure 
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L = 2.5 nH. 



The scattering matrix of an impedance Z connected in series with the trans- 
mission line is : 



S = 



■^11 

^21 



•^12 

*^22 



2+2 2+z 

2 z 

2+z 2+z 



(4.21) 



where 



2o ’ 



and Zo is the characteristic impedance (50 Ohm) of the connected line. In this case, the 
TOUCHSTONE CAD program is useful in calculating the scattering parameters. 

These theoretical outputs are shown in Appendix C. Also the magnitude plots 
of the calculated scattering parameters of this device are shown in Fig. 15. 
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Figure !5. Tlieoretlcal Scattering Parameters of 51 Ohm Resistor 
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2. Device Embedment 

The S-paramcters of both fixtures plus device are measured by the network an- 
alyzer at planes AA and BB” as shown in Fig. 16. Typical measured magnitude data is 
shown in Fig. 17 with data tables in Appendix C. 



INPUT HALF 
OF FIXTURE 



DEVICE 



OUTPUT HALF 
OF FIXTURE 





Netv\/ork analyzer 
measurement plane 




Figure 16. Signal Flow graph for the Measurement of the (Fixture + Device) S- 
parameters 
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Figure 17. The Magnitude Plot of Both Microstrip Fixture and 51 Ohm Resistor 
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3. De-embedding the Device 

Now we are ready for calculation of the S-parameters of the device. The S- 
parameters of fixture A, fixture B and both fixtures plus device constitute the necessary 
data. Once the microstrip circuit board is prepared, then the network analyzer provides 
the three different input reflection coefficients F ,„2 and F^j with respect to three dif- 
ferent calibration loads. Also computations are carried out to find the load reflection 
coefficients of these three different loads, F^,, F^j and F^j. After following the same 
procedure for both fixture A and fixture B, we can get the respective S-parameters. 
From the matrix equation (Eq. 2.17), the three unknowns S,,, S 22 and A are solved. For 
de-embedding the device itself, we need the additional S-parameters for the device plus 
fixture. These three sets of S-parameters are utilized to solve the de-embedding problem. 

Equations 4.11 to 4.14 are now ready for use. Through these equations, the 
de-embedded outputs are obtained, as shown in Fig. 18 (the table of de-embedded S- 
parameters is shown in Appendix C). In addition to the plot of the de-embedded mag- 
nitudes, the comparison of the scattering parameters Sn, S ,2 and S 22 of the de-embedded 
output with the parameters of the embedded resistor and the theoretical values are also 
available in Fig. 19 to Fig. 21 (data tables in Appendix C). 
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Figure 18. Magnitude Plot of tlie De-enibcdded Parameters of 51 Ohm Resister 
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Figure 19. Comparison of the Sll of the 51 Ohm Resistor «ith the Theoretical 
Values 
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Figure 20. Comparison of the S12 of the 51 Ohm Resistor «ith the Theoretical 
Values 
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Figure 21. Comparison of the S22 of tlie 51 Ohm Resistor with the Iheoretical 
Values 
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V. FURTHER DE > EMBEDDING METHODS 



A. MEASUREMENT OF DEVICE PARAMETERS USING A SYMMETRIC 
FIXTURE. 



1. Background 

As discussed earlier, if the device is constructed without a coaxial or waveguide 
interface, de-embedding becomes more difficult. In some cases, it is possible to proceed 
using an auxiliary standard whose properties are calculated or measured at low fre- 
quencies, but, in other cases, no suitable auxiliaiy’ standard is available. In either case, 
if the fixture is symmetric, a de-embedding procedure that does not require the use of 
an auxiliary standard can be developed [Ref 5]. 

The concept for symmetric de-embedding is illustrated in Fig. 22. Fig. 22-a 
shows a device inserted in a symmetric test fixture. We assume that the device has 
complementary interfaces on the input and output sides so that when the device is re- 
moved from the fixture, the two halves of the fixture can be connected, as shown in Fig. 
22-b. The scattering parameters of each half of the test fixture can be extracted from the 
scattering parameters of the whole test fixture connected, as illustrated in Fig. 22-b. 
Then the device is inserted into the test fixture, as illustrated in Fig. 22-a, and the scat- 
tering parameters are again measured. 

The method for implementing the concept outlined above is as follows. If we 
have a fixture whose T-parameters in a 50 ohm system are CTf], then if CZ/Tfl re- 
presents the square root of CTF], 



Now suppose we want to measure the T-parameters of a device in a 50 ohm system 
[rS]. If we separate the fixture into two symmetric parts and inserted the device to be 
measured, then for the combination we would measure a set of T-parameters CTFS] 
given by: 






(5.1) 



CTFS: = lHTr\lTS~\[.HTF:\. 



(5.2) 



To find CTS], we would simply invert above equation to get 

CrS] = IHTFT^ITFS:\IHTFT^. 



(5.3) 
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( TF 1 (b) 

Figure 22. Device Embedded Between Symmetric Fixture and Test Fbcture 



This procedure de-embeds the device parameters CFS] from the measured parameters 

ZTFS'] and CTF]. The steps required to find CFS] in a symmetric fixture are as follows: 

1. Measure the scattering parameters of both fixture plus device, C5FS], using the 
network analyzer, 



ZSFSJ = 



SFS,, 

SFS21 



sfs ,2 

SFS22 



(5.4) 



then transform to the T-parameters: 



crFS] = 



SFS2X 

-SFS22 

SFS21 



SFSu 

SFS2: 

1_ 

SFS2, 



(5.5) 



where iV/, = SFS,,SFSj 2 + SFS 12 SFS 21 . 

2. Measure the scattering parameters of both fixtures CSF]. In practice, both fixtures 
are measured with the network analyzer individually: 
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iHSn = [//SFl] = \:HSF22 = 

Then, transforming to the T-parameters: 

-A/, //Sf,, 



7757,, 


7757,2 


77572, 


775722 



\: htf 2 = 



HSF- 



21 



HSF 



21 



-!^SFn 

HSF 



21 



HSF 



21 



where .V/^ = HSF^HSF^^ + HSFJiSF.y. 
3. Invert £HTF2: 



IHTF 2~^ = 



HTF, 



11 



M , 

-HTF, 



12 



- HTF^x 

A/3 

HTF 



11 









where M, = HTF,, HTF,, - HTF„HTF„. 
4. Find device T-parameters CFS]: 



[7-5: = 



5. Transform FT'S] to S-parameters 
•^11 

Sj , S22 



■7777,, 


7777,2' 


-1 


'TK,, 


775,2' 




■7777,, 


7777,2' 


.77772, 


7/7722_ 




_rfs„ 


77522_ 




.77772, 


T77722_ 



CSS] = 
where 



^ 7-5,2 

^11 - 7-5 
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5i 2 — T5,j — 



75 , 2752 , 



75 



22 



■^21 ~ ' 



522 ~ “ 



1 

7522 
75 



75. 



21 

22 



(5.6) 

(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 
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2. De- Embedding of a Filter Unit 

As a test demonstration of the symmetrical de-embedment procedure, an 8436A 
bandpass filter (8-12 GHz) is used as the test device. All components of the network 
had type APC-7 test port connectors. The outputs [SFS] and ZHSF2 of the network 
analyzer are shown in Appendix D. To get the de-embedded outputs, the Fortran pro- 
gram shown in Appendix F is used. The magnitude plot of the filter embedded, the filter 
itself, and the de-embedded filter are shown in Figs. 23 to 25. Data Tables are shown 
in .Appendix D. The comparison of the scattering parameters 5,„ and Sj, of the filter 
itself, and embedded filter are shown in Figs. 26 to 28. Tabulated data is shown in Ap- 
pendix D. 



40 



dB 

Loss 




Figure 23. Tlie S-paranieter Magnitude Plot of Both Fixture and Filter 
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Figure 24. The S-paranieter Magnitude Plot of Filter Without Fixtures 
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Figure 25. Tlie S-paranieter Magnitude Plot of De-embedded Parameters of Filter 
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Figure 26. Comparison of the Sll of De-enibedded Filter Mith Filter \\ithout Fix- 
tures 
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Figure 27. Comparison of the S 12 of De-einbedded Filter «ith Filter Without Fix- 
tures 
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B. TSD (THROUGH, SHORT, DELAY) METHOD 

This procedure, developed by R. Speciale [Refs. 6, 7 ], is intended for system calibration 
and systematic error removal in automated S-parameter measurement. This procedure 
uses only three non-critical reference standards, called the "Through", the "Short" and 
the "Delay", and leads to a set of equations for which a fast, non-iterative, e.xplicit sol- 
ution is known [Ref 6]. 

1. The Basic Concept 

This procedure uses only three mechanically simple reference standards, known 
respectively as the "Through", the "Short", and the "Delay". 

The "Through" is defined as the direct mutual connection of fixtures A and B. 
as shown in Fig. 29, when the uncalibrated S-parameter data for the device are acquired. 
The S-parameter data for "Through" are therefore acquired with the measurement lines 
to the device directly interconnected and the device removed from the signal path. 

The "Short" is defined as a reference object presenting perfect shorts at the im- 
mediate interfaces where the unknown device is connected for measurement. 

Finally, the "Delay" is defined as any arbitrary length of 50 ohm line inserted 
between the interfaces where the unknown device is connected for measurement. Sum- 
marizing. for the "Through" and the "Delay", the signal path is continuous between ports 
1 and 2 of the S-parameter unit, through the couplers used for connecting the object. 
•A direct connection represents the "Through" and a finite arbitrary length line for the 
"Delay", in place of the unknown object. On the contrar>‘ the signal path is blocked for 
the "Short" by perfect shorts appearing at the interfaces where the unknown device is 
normally connected. 

Three S-parameter matrices II-Sj] and are acquired at each desired 

frequency during system calibration. These matrices describe the "Through", the "Short" 
and the "Delay" respectively. Because of the device imperfections, these matrices do not 
represent theoretically perfect "Through", "Short" and "Delay" elements, which are as- 
sumed to exist in the actual measurement. 

Now consider the error boxes A and B. Error box A is assumed to be connected 
on the side of port 1, while error box B is assumed to be on the side of port 2, as shown 
in Fig. 29. 
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NETWORK - ANALYZER 




Figure 29. Block Diagram of the TSD Calibration Procedure 



From Fig. 29, the error boxes A and B are always assumed to be inserted at the 
two ends of the measurement 50 ohm lines connecting the device to ports I and 2 of the 
S-parameter unit. An algorithm has been derived wh -h performs an explicit, non- 
iterative calculation of the S-parameter matrices ^nd which characterize the 

virtual error bo.xes A and B. 

2. The Equations and their Explicit solution 

As discussed earlier, the T-matrices can be multiplied simply with cascaded ele- 
ments. Therefore, consider Fig. 30, where the T-matrix for the connection shown is: 

= ITsl. (5.14) 

From the above, we get: 

cro] = cr^r' cTsr' 
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(5.15) 
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Figure 30. Device Embedded Setup 



where CFo] is the corrected T-matrix and Cr^o^Il is the uncalibrated T-matrix for the 
device. 

The T-parameters of fixtures plus device, Jire given by tlie network 

analyzer. But the problem is to obtain the T-parameters of each fixture, CT<II and 
CTb]. To do this, we consider the error boxes A and B as shown in Fig. 31. 

From a formal point of view, both the "Through" and the "Delay" can be con- 
sidered as different, finite lengths of 50 ohm lines, cascaded between the error bo.xes A 
and B. With CT^iD and the T-matrix derived from the S-parameter data obtained 

in these two calibrations, and CT,(I1 and CTbU, the unknown T-parameters of the error 
boxes A and B, we can write: 

= (5-16) 

For the 50 ohm lines L, and Li we have; 






,-pl 



0 



0 






(5. IS) 



where +J^^. 
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Figure 31. Error Boxes of A and B with Delay 



cr^2: = 





(5.19) 



where Pi = r^+ jOi- 

Equations (5.16) and (5.17) together with the two equations for the "Short" 
characterize boxes A and B. Also, consider the situation for the "Short", as shown in 
Fig. 32, we have 



S 



115 ~ 






I + •^22.4 



(5.20) 



5 



115 ~ 



^28 



^228 ~ 1 ~ 

1 + •^ns 



(5.21) 



where S,,j and S 22 S are S„ and S 22 of the "Short", and S„^, S,,^ are S„ of box A and box 
B (det means matrix determinant), where i = 1 . 2. 

From Eqs. (5.16) and (5.17), we get: 



n 2.4 

■^1 1.4 ~ T 

‘22A 



Hi 

^1 



(5.22) 



50 




(ict — *^22/1 



A 

H, 



S-y-tp — ' 






^22B - r ~ i<r 
h2B 



-K, 



det Sq — S] 13 ^ 



Now substitute Eqs. (5.22) and (5.22a) into (5.20) to obtain: 



^22A - 



H, 



' ~ -^ii: 



H, 



and also, 



H, 



^22A 



(5.22a) 

(5.23) 

(5.23a) 



(5.24) 



(5.25) 
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Substituting Eqs. (5.23) and (5.23a) into (5.21), 



A-a 

r ~^22S 






^22S+ 



(5.26) 



~ •^'uB-^22S ~ 



^22B + 



A', 



(5.27) 



>115 



where 



^ 2 L - P J 






1 +■ 



^22 -An 
2 sinh p 



] 



(5.28) 



(5.29) 






12 



2 sinh p 



A, = - 



A, 



12 



2 sinh p 



H, = - 



Hix 



7 <;1 



sinh p 



A,= 



A 



21 



4 -) 



2 sinh p 



(5.30) 



(5.31) 



(5.32) 



(5.33) 



Here the quantities and A)^ are the elements of the matrices CZZ3 and CA], defined 
bv 



c//: = C7'o,3 UdiT' 

CA: = cro2r' 



(5.34) 

(5.35) 



and p is defined by; 



p cosh ^ 2 (^11 + ^22)^ ~ cosh ^ ^ (^11 "h ^22)^ 



(5.36) 



Once the S-parameters of the error boxes A and B arc known as given by Eqs. (5.22) 
through (5.27), use the general conversion formulas, 





Tn 


1 


■-dots 


•^11 


7^1 


1 

t-r* 


^ 2 \ 


•^22 


1 _ 



(5.37) 



where det 5 = — -Sjj.Sjj. 

As discussed before, the following calculation to de-embed the device matrix, 

cr^] = cr^r‘ cr^r*. 

can be carried out using Eqs. (5.22) through (5.27). 

3. Discussion 

This calibration procedure uses only three reference standards, which is the 
minimum necessary' for characterizing the error boxes A and B. By using two delays of 
diOerent length, for a total of 4 standards, a check of intrinsic consistency may be ob- 
tained. The standards used are simple enough to be inserted at the location of the ob- 
ject, even if this is an active device contained in some kind of measurement fixture. In 
this way, all system errors are removed, including those introduced by the fixture. 

The position of the reference plane is immaterial, as long as it is kept constant 
during calibration and measurement. Also, the given solution of the relevant network 
equations is an exact, explicit solution. Finally, the procedure is essentially independent 
of the specific system configuration. 
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VI. CONCLUSION 



Various sources have been cited for the extraction of the scattering parameters of a 
device by removal of the error and contributions of the fixturing setup [Refs. 8,9.10,111. 
In Chapters II and III, the de-embedding procedure using microstrip line within the 
given conditions was discussed. This project included various aspects such as the 
microstrip line fixturing setup and the construction and handling of small experimental 
boards. But in spite of the limitations of the early-model Network analyzer HP 8409B 
which w'as available for this work, the best possible results were obtained using repeated 
measurements on the same experimental board to obtain mean values. 

To extract the scattering parameters of the device 51 ohm resistor, the nine basic 
parameters as discussed earlier in Chapter IV (Eqs. 4-11 to 4-19) were used. These 
equations required the scattering parameters of both fixtures since they are not sym- 
metrical. That is, individual measurements and calculations were required for both fix- 
tures. The matrix equation (Eq. 2-17) involved the three unknown scattering parameters 
S,„ $22 and A = (SnSjj — -SuSsi) for fixture A and B. That meant each fixture required 
three equations for the matrix solution. Three different measurements and calculations 
of microstrip line were used to set up the above matrix entries, i.e, three different meas- 
ured reflection coefficients from the network analyzer with respect to the three different 
microstrip loads, respectively r„„ and r,„j and three different calculated reflection 
coefficients at the terminals of the three different loads fn, and fi ,3 were required. 
The computed value for S„ and S 22 are combined with a knowledge of A which can be 
used to determine the SijSj, product. Since S, 2 , and Sj, for the fixtures always appear in 
product form, the magnitude of these terms can be calculated as the square root of the 
product term in the expression for A . 

The de-embedded scattering parameters of the device (51 ohm resistor) were found 
through experiment to be close to the theoretically calculated values, shown in magni- 
tude plots in Figs. 19, 20, and 21 and as tabulated in Table 4. As may be seen in Figs. 
19 to 21, the theoretically calculated scattering coefficients of the 51 -Ohm resistor device 
showed only a slight continuous variation over the entire frequency range of interest. 
The measured values of the embedded device showed, however, rapid irregular fluctu- 
ations with frequency, reflecting the probable occurrence of internal resonances within 
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parasitic discontinuities of the coupling circuitr>\ It is the task of the de-embedding 
procedure to remove these effects from the measurement. The curves numbered (3) in 
the above Figures show the de-embedded measurements, in comparison with the the- 
oretical curves. As may be seen, the de-embedded results followed the e.xpected values 
closely, showing fluctuations which are assumed to be due to measurement error, and 
network-analyzer calibration errors. 



Table 4. THE COMPARISON OF THE OUTPUT FOR THE 51 
OHM RESISTOR IN DB UNITS 



Freq.(MHz) 


Theoretical 


De- 

embedded 


Embedded 


2000.0 


6.11 


6.35 


11.14 


2300.0 


5.54 


5.12 


10.18 


2600.0 


5.02 


4.22 


7.13 


2900.0 


4.55 


3.87 


6.21 


3200.0 


4.13 


3.19 


7.33 


3500.0 


3.76 


3.16 


11.29 


3S00.0 


3.42 


3.91 


11.59 


4100.0 


3.12 


2.66 


6.27 


4400.0 


2.86 


4.11 


6.07 


4700.0 


2.62 


2.25 


9.91 


5000.0 


2.41 


2.11 


11.95 



5,2 and Sji terms for fixture A and B always appears in product form, considering 
their phase is not important. As long as the pliase angles of S 21 for fixture A and Sj, for 
fixture B are approximately the same, one can choose either of the phase angles. Some 
confusion may arise, if the scattering parameters of both fixture A and B are specified 
independently of the final device (51 ohm resistor) parameters. In this case, the ap- 
proximate phase shift through one fixture can be computed, or determined exper- 



imentally, and used to choose the proper phase of Sj, for fctture A and Sj, for fixture 
B. 



In Chapter V, the de-embedding procedure with symmetrical fixtures was carried 
out, using as the device a commercial bandpass filter (8-12 GHz). Standard hardware 
couplers were added to represent the fixtures. Since the fixtures were symmetrical, the 
de-embedding procedure was relatively straightforward. As may be seen in Figs. 26 to 
28. the de-embedded scattering parameters of the filter were very close to the expectation 
values, i.e. of the scattering parameters of filter without fixtures. It may also be noted 
that the measured S-parameters of the embedded filter (fixtures plus filter) were very 
similar to those of the filter alone, except of course for phase; indicative of well-matched 
coupling hardware. Therefore, in this task, the de-embedded results were found to be in 
good agreement with the expected values, as shown in Figs. 26 through 28. 
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APPENDIX A. THE OUTPUTS OF THE NETWORK ANALYZER AND 
THEIR REFLECTION COEFFICIENTS FOR FIXTURE A AND B 




ur-'s ■ nnr 



Figure 33. Measurement of the Load 1 (Fixture A) 



Table 5. THE MEASUREMENT OF THE LOAD KFIXTURE A) 



Freq.l MHz) 


dll (lo?«;) 


.Ansletdeg. ) 


r„ 


2000. U 


3.07 


T-> 7 


0.65-^ jo. 27 


2150.O 


3.82 


-63.1 


0.29 - jo.57 


2300.0 


3.79 


-153.9 


- 0.58 - jo. 2 8 


24.^0.0 


3.46 


126.4 


.0.4(1 •+-](). 54 


2600.0 


3.20 


60.5 


0.34 -(- ii i.60 


2750.0 


3.33 


-1.7 


0.6S - jo. 02 


29O0.O 


4.04 


-69.9 


0.22 - jO.59 


?05n,o 


5.39 


-152.8 


-0.4S - ]0.25 


3200.0 


6.17 


107.6 


- 0 . 15-^jo.47 


3350.0 


5.32 


17.1 


0.52 + jo. 16 


3500.0 


4.65 


-51.0 


0.37 - jo.45 


3650.0 


3.97 


-111.9 


-0.24 - jO.59 


3800.0 


4.89 


-175.9 


-0.57 - ]0.04 


3950.0 


7.73 


104.5 


-0. 1 0 -4- jo.4o 


4100.0 


7.36 


-17.1 


0.41 - j(».13 


4250.0 


6.70 


-121.4 


-0.24 - jo.39 


4400.0 


4.59 


1 70.4 


-0.58 jo. 10 


4550.0 


4.12 


1 14.2 


- 0 . 26 -^ jo .57 


4700.0 


4.88 


58.9 


0.29-t- jo.49 


4S50.0 


7.91 


-9.1 


0.40 - j 0.06 


5000.0 


8.42 


-137.9 


-0.28 - jO.25 
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Figure 34. Measurement of the Load 2 (Fixture A) 



Table 6. THE MEASUREMENT OF THE LOAD 2 (FIXTURE A> 



Frcq.fMI Iz) 


dB (los<:) 


Angle! deg.) 


r,„ 


2<)<)0.i» 


2.28 


98.2 


- 0. 1 1 4“ jo. 76 


21.50.0 


2.89 


29.5 


0.62- |0. 3 > 


2300.0 


3.77 


-49.2 


0.42 - )0.49 


2450.0 


4.12 


-138.0 


-0,46 - jU.42 


2600.0 


3.75 


136.6 


-().47 'h ]0.45 


2750.0 


3.31 


65.1 


0.2^-^j0.62 


2900.O 


3.40 


4.9 


0,67 -r j().06 


3050.0 


3.60 


-54.2 


U.39 - ]0.54 


3200.O 


4.55 


-119.0 


-0.29 - jo. 52 


3350.0 


5.78 


159.8 


-0.48 -]0. IS 


3500.0 


6.39 


61.8 


0.23 4- ]0.42 


3650.0 


5.76 


-25.5 


0.47 - jo.22 


3800.0 


5.70 


-89.0 


0.01 - jO.52 


3950.0 


4.69 


-134.2 


-0.41 - jo.42 


4100.0 


3.91 


173.7 


-0.63 + j0.07 


4250.0 


4.48 


117.4 


-0.27^)0.53 


4400.0 


6.84 


47.9 


0.31 +j0.34 


4550.0 


S.48 


-61.4 


0.18 - )0.33 


4700.0 


7.57 


-174.8 


-0.42 - jo.04 


4850.0 


4.92 


115.1 


-0.24^)0.51 


5000.0 


3.99 


62.0 


0.30-!- ju.56 
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Figure 35. Measurement of the Load 3 (Fixture A) 



Table 7. THE MEASUREMENT OF THE LOAD 3 (FIXTURE A) 



r-rcq.( MHz) 


(JB (loss) 


.*\nele(dcg.) 


r„ 


2UUO.O 


2.14 


175.1 


-0.7S-^jO.07 


2150.0 


2.21 


114.5 


-0.32^|O.7| 


23OU.0 


2.57 


55.4 


U.42-J0.6I 


2450.0 


3.25 


0.1 


0.69 - jo.06 


2600.0 


3.86 


■74.3 


0. 1 7 - ]0.62 


2750.0 


4.44 


-156.0 


-0.55 - jo. 24 


2900.0 


4.42 


119.8 


-0.3o- 10.52 


30.5O.O 


3.9S 


46.6 


0.43 ^46 


3200.0 


3.89 


-14.2 


0.62 - 


3350.0 


3.S9 


-66.5 


0.25 - jo. 59 


3500.0 


4.15 


-120.7 


-0.32 - jo.53 


3650.0 


5.02 


175.9 


-0.56-HjO.O4 


3S00.0 


8.22 


99.0 


-0.06-HjU.38 


3950.O 


8.95 


9.5 


O.35-^j0.Ob 


4 1 oO.O 


6.37 


-79.5 


0.09 - jo.47 


4250.0 


4.35 


-147.2 


-0.51 - jO.33 


4400.0 


3.49 


157.8 


-0.62-j0.25 


4550.0 


3.37 


108.4 


-0.21 jO.64 


4700.0 


3.90 


58.8 


0.33 -H jO.55 


4850.0 


5.55 


1.4 


0.53 -H jo. 01 


5000.0 


6.59 


-82.6 


0.06 - jO.46 
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Figure 36. Measurement of the Load I (Fixture B) 



Table <S. THE MEASUREMENT OF THE LOAD I (FIXTURE B) 



F'req.(MlIz) 


dl( (loss) 


■Angleidee.) 




2000.0 


1.97 


-ir.s 


-0.37 - ]0.71 


2150.0 


1.S9 


-1"6.2 


-O.SO - ]0.()5 


2300.0 


2.00 


130.2 


-0.51 ]0.61 


2450.(J 


2.27 


81.9 


0.1 1 + j0.7() 


26O0.0 


2.40 


34. S 


0.62 + )0.43 


27.50.0 


2.62 


-13.9 


0.72 - jo. IS 


29< ) 0 .o 


3.04 


-6". 3 


0.27 - jo.65 


3(.).50.0 


3.S7 


-130.3 


-0.41 - jo.49 


320O.O 


5.05 


151.6 


-0.49 4-jo.27 


3350.0 


5.41 


64.0 


0.24-r jo.48 


3500.0 


4.91 


-11.^ 


0.56 - jo. 12 


3650.0 


3.87 


-67.6 


0.24 - jo.5v 


3S00.0 


2.90 


-118.9 


-0.35 - jO.63 


3950.0 


2.55 


- 1 66.3 


-0.72 - jo. 18 


4100.0 


2.68 


147.0 


-0.62 + jo. 40 


4250.0 


3.33 


96.3 


-0.07 -h j().6S 


4400.0 


4.95 


36.4 


0.46 4- jo. 34 


4550.0 


7.33 


-•^3.-1 


0.31 - j<».3(» 


47oo.O 


7.51 


-139.- 


-0.32 - jo.27 


4S5O.0 


5.35 


143.6 


-0.43 + j0.32 


5000.0 


3.67 


87.9 


0.02 + jo.65 
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Figure 37. Measurement of the Load 2 (Fixture B) 



Table THE MEASUREMENT OF THE LOAD 2 (FIXTURE Bl 



Freq.t MHzi 


clB (loss) 


.•\ncie(dec. ) 


r,„ 


2t*(»0.0 


1 J 
1 J 


-26.9 


0.68 - )0.35 


2150.0 


1 J 


-91.9 


-0.02 - (0.75 


230O.O 


2.44 


-155.5 


-0.60 - ji ».3 1 


245O.0 


2.51 


147.3 


-o.63^)0.4(i 


2600.O 


2.41 


96.6 


-0.0M-jo.75 


2750. 0 


2.25 


50. 1 


0.50-t- jo.5'^ 


290O.O 


2.35 


6.9 


0.76-Hjn.()y 


3i.)5o.o 


2.2S 


-35.9 


0.62 - ju.45 


320O.O 


2.51 


-7S.5 


0.15 - |0.73 


3350.0 


2.69 


-123.9 


-0.41 - j0.6l 


350U.O 


3.25 


-175.1 


-0.69 - ]0.0(' 


3650.0 


4.27 


126.2 


-0.36-<-i().49 


3SOO.0 


5.24 


56.1 


0.3 1 ]0.45 


3950.0 


5.38 


-20.6 


0.50 - jO.l') 


4100.0 


4.49 


-90.2 


0.(«,) - j0.6<» 


4250.0 


3.49 


-147.9 


-0.57 - j0.3o 


4400.0 


2.76 


163.0 


-0.7o-jO.21 


4550.0 


2.43 


11S.9 


-0.37 jo.oo 


4700.0 


2.36 


77.5 


0.l6^jO.74 


4850.0 


2.55 


36.2 


0.60-t-)0.44 


5000.0 


3.32 


-S.6 


0.67 -jO. 10 
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Figure 38. Measurement of the Load 3 (Fixture B) 



Table 10. THE MEASUREMENT OF THE LOAD 3 (FIXTURE B) 



Frcq.(MHz) 


dB ( loss) 


.-\nclc(deg.) 


r. 


2< KJU.U 


1.41 


78.5 


0.17-j.),S3 


2150.0 


1.71 


28.1 


0.72^)0.39 


230<.i.<> 


2.16 


-20.2 


()."() . |ti,34 


2450.0 


2.4S 


-84.3 


0.0“ - 


26OO.0 


2.69 


-145.9 


-0.61 - io.41 


2750.0 


2.85 


152.4 


-0.64-- jO.33 


2900.0 


2. 85 


94.8 


-0.06-^-|0.72 


.'O50.O 


2.75 


43.2 




320O.0 


77 


-2.^ 


0.7' - j0.o4 


3350.0 


2.68 


-44.3 


0.53 - ]<k51 


35(;o.o 


2,68 


-82.4 


0. In - 10.73 


3650.0 


1.93 


-121.1 


-0.41 - jo. 69 


3SO0.0 


1.91 


-162.1 


-0.76 - io.25 


3950.0 


2.02 


157.0 


-0.73 + )0.31 


410O.0 


2.42 


113.2 


-().30 4- jo.7() 


4250.0 


3.21 


64.6 


0.30-^ jo.62 


4400.0 


4.49 


7.5 


0.54-^ |(M)S 


4550.0 


5.70 


-60.8 


0.2.5 - jo. 45 


4700.0 


5.80 


- 1 .>4.9 


-0.3(' - jo. 3 6 


4850.0 


4.89 


158.8 


\ 

i 


5000.0 


3.72 


105.1 


-0. 17 + jO.63 
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APPENDIX B. THE S-PARAMETERS OF THE FIXTURES 



Table 1 1. THE S - PARAMETERS OF FIXTURE A 



Freq.(MHz) 


Au 


An 


A^2 


2000.0 


-0.130 j0.0S2 


0.837 - jO.143 


0.132 j0.029 


2150.0 


-0.064 + jO.153 


0.668 - jO.499 


0.139 - j0.019 


2300.0 


0.026 -i- j0.165 


0.355 - jO.738 


0.130 - jO.OSO 


2450.0 


0.099 + j0.125 


0.021 -i- jO.802 


0.100 - j0.147 


2600.0 


0.141 -(- j0.069 


0.380 + jO.701 


0.034 - j0.17S 


2750.0 


0.151 - j0.007 


0.653 4- jO.443 


-0.018 - jO.lSO 


2900.0 


0.091 - j0.063 


0.777 -1- j0.0S7 


-0.032 - j0.128 


3050.0 


0.048 " j0.055 


0.717 • jO.287 


-0.003 - j0.074 


3200.0 


0.042 - jO.026 


0.485 - jO.576 


0.060 - j0.037 


3350.0 


0.098 - jO.025 


0.151 - jO.720 


0.180 - jO.031 


3500.0 


0.137 - jO.113 


0.199 + jO.679 


0.288 - jO.088 


3650.0 


0.054 - jO.265 


0.470 4- jO.482 


0.330 - j0.184 


3S00.0 


-0.148 - j0.261 


0.579 4- j0.214 


0.327 - jO.254 


3950.0 


-0.258 - jO.163 


0.625 - j0.036 


0.321 - jO.239 


4100.0 


-0.334 -i- j0.061 


0.581 - jO.305 


0.299 - jO.322 


4250.0 


-0.256 + jO.288 


0.318 - jO.506 


0.442 - jO.181 


4400.0 


0.174 4- j0.317 


0.008 - jO.477 


0.713 - j0.044 


4550.0 


-0.502 - jO.872 


0.225 - jl.580 


1.743 - jl.374 


4700.0 


-0.064 - j0.531 


0.670 - j0.309 


0.034 - j0.761 


4850.0 


0.183 + jO.229 


0.651 + jO.057 


0.137 - j0.304 


5000.0 


0.062 + j0.136 


0.646 - jO.289 


0.219 - jO.013 
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Table 12. THE S - PARAMETERS OF FIXTURE B 



Frcq.(MHz) 


Bn 


Bn 


Bn 


2000.0 


-0.060 + jO.107 


0.342 - jO.834 


0.036 - j0.076 


2150.0 


-0.007 + j0.116 


0.057 - jO.892 


0.008 - jO.081 


2300.0 


O.027 + jO.lOO 


0.230 + j0.850 


-0.032 - j0.085 


2450.0 


0.049 + jO.091 


0.475 + jO.723 


-0.083 - j.082 


2600.0 


0.080 + jO.084 


0.669 “i" jO.533 


-0.115 - jO.075 


2750.0 


0.124 + jO.063 


0.792 - j0.291 


-0.125 - jO.062 


2900.0 


0.152 + j0.007 


0.831 + jO.028 


-0.136 - jO.029 


3050.0 


0.158 - jO.073 


0.788 - jO.229 


-0.125 - jO.007 


3200.0 


0.106 - j0.143 


0.669 - jO.451 


-0.110 + jO.016 


3350.0 


0.026 - jO.168 


0.494 - j0.630 


-0.064 + j0.036 


3500.0 


-0.035 - j0.149 


0.268 - jO.741 


0.016 + j0.059 


3650.0 


-0.092 - jO.112 


0.033 - j0.790 


0.108 + j0.019 


3800.0 


-0.076 - jO.072 


0.230 + jO.762 


0.165 + j0.018 


3950.0 


-0.069 - jO.097 


0.468 + jO.633 


0.261 + j0.0I9 


4100.0 


-0.137 - j0.157 


0.636 + jO.412 


0.389 + jO.014 


4250.0 


-0.331 - j0.130 


0.681 + jO.135 


0.543 + jO.018 


4400.0 


-0.572 + jO.161 


0.478 - jO.116 


0.808 + jO.024 


4550.0 


-0.114 + jl.l7l 


1.062 4- jO.589 


1.625 - jO.527 


4700.0 


0.230 - jO.171 


0.577 - jO.759 


0.163 - jO.968 


4850.0 


-0.065 + j0.166 


0.306 - jO.676 


0.219 - jO.546 


5000.0 


0.124 4- jO.326 


0.091 - jO.685 


0.291 - jO.446 
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APPENDIX C. THE S-PARAiMETERS OF THE EMBEDDED AND 
DE-EMBEDDED 51 OHM RESISTOR 



Table 13. THE THEORETICAL S - PARAMETERS OF THE 51 
OHM RESISTOR IN DB UNITS. 



Freq.(MHz> 


•^11 


•S',: 


*^21 


5:: 


2000.0 


6.11 


4.27 


4.27 


6.11 


2150.0 


5.82 


4.37 


4.37 


5.82 


2300.0 


5.54 


4.47 


4.47 


5.54 


2450.0 


5.28 


4.58 


4.58 


5.28 


2600.0 


5.02 


4.69 


4.69 


5.02 


2750.0 


4.78 


4.81 


4.81 


4.78 


2900.0 


4.55 


4.93 


4.93 


4.55 


3050.0 


4.34 


5.05 


5.05 


4.34 


3200.0 


4.13 


5.17 


5.17 


4.13 


3350.0 


3.94 


5.30 


5.30 


3.94 


3500.0 


3.76 


5.43 


5.43 


3.76 


3650.0 


3.58 


5.56 


5.56 


3.58 


3S00.0 


3.42 


5.69 


5.69 


3.42 


3950.0 


3.27 


5.82 


5.82 


3.27 


4100.0 


3.12 


5.95 


5.95 


3.12 


4250.0 


2.99 


6.09 


6.09 


2.99 


4400.0 


2.86 


6.22 


6.22 


2.86 


4550.0 


2.73 


6.36 


6.36 


2.73 


4700.0 


2.62 


6.49 


6.49 


2.62 


4S50.0 


2.51 


6.63 


6.63 


2.51 


5000.0 


2.41 


()."(> 


6.76 


2.41 



65 



Table 14, THE S - PARAMETERS OF BOTH FIXTURE AND 51 
OHM RESISTOR 



Freq .( MHz ) 








A /:. 


2000.0 


11.14 


4.40 


4.60 


6.45 


2150.0 


11.38 


4.64 


4.71 


6.27 


2300.0 


10.18 


4.90 


4.80 


6.05 


2450.0 


8.47 


5.36 


5.35 


5.86 


2600.0 


7.13 


5.87 


5.84 


5.51 


2750.0 


6.35 


6.30 


6.25 


4.99 


2900.0 


6.21 


6.99 


6.99 


4.68 


3050.0 


6.51 


7.54 


7.49 


4.48 


3200.0 


7.33 


7.32 


7.25 


4.43 


3350.0 


8.92 


7.04 


6.65 


4.94 


3500.0 


11.29 


7.10 


5.81 


5.75 


3650.0 


12.01 


6.46 


6.31 


6.18 


3800.0 


11.59 


7.30 


7.25 


6.91 


3950.0 


9.15 


6.45 


6.90 


7.83 


4100.0 


6.27 


6.35 


6.56 


9.07 


4250.0 


5.53 


6.00 


6.25 


1.78 


4400.0 


6.07 


5.61 


5.85 


12.72 


4550.0 


7.54 


5.30 


5.50 


12.53 


4700.0 


9.91 


5.34 


5.39 


10.00 


4850.0 


10.46 


5.70 


5.66 


7.82 


5000.0 


11.95 


6.10 


6.05 


6.53 
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Table 15. THE DE - EMBEDDED OUTPUT OF 51 OHM RE- 
SISTOR DEVICE UNDER TEST. 



Frcq.(MHz) 


‘^11 


5.: 




5:: 


2U00.0 


6.35 


3.62 


3.82 


4.92 


2150.0 


5.72 


3.97 


3.64 


5.84 


2300.0 


5.12 


3.97 


3.97 


4.74 


2450.0 


4.57 


3.18 


3.17 


3.75 


2600.0 


4.22 


3.39 


3.36 


3.04 


2750.0 


3.61 


3.66 


3.61 


4.74 


2900.0 


3.87 


4.27 


4.27 


3.77 


3050.0 


3.64 


4.91 


4.86 


2.91 


3200.0 


3.19 


4.51 


4.44 


2.67 


3350.0 


3.03 


4.16 


3.77 


2.41 


3500.0 


3.16 


4.72 


4.43 


3.79 


3650.0 


3.62 


4.65 


4.50 


4.45 


3S00.0 


3.91 


4.56 


4.51 


3.72 


3950.0 


3.35 


5.38 


4.83 


2.73 


4100.0 


2.66 


5.06 


5.27 


4.66 


4250.0 


3.45 


5.71 


5.46 


5.31 


4400.0 


4.11 


5.83 


6.07 


2.52 


4550.0 


3.17 


6.36 


6.56 


2.07 


4700.0 


2.25 


6.73 


6.68 


3.34 


4850.0 


2.08 


6.73 


6.69 


2.95 


5000.0 


2.11 


6.96 


6.91 


2.02 
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APPENDIX D. THE S-PARAMETERS OF THE FIXTURE, FILTER 
EMBEDDED, FILTER WITHOUT FIXTURES AND FILTER 

DE-EMBEDDED 



Table 16 . THE S - PARAMETERS OF THE FIXTURE 



Freq.(MHz) 


■s.. 






6000.0 


-0.08 + j0.17 


0.81 + jO.55 


-0.09 4- j0.13 


6320.0 


-0.09 + j0.04 


0.76 - j0.61 


-0.17 + jO.lO 


6640.0 


-0.06 + j0.04 


-0.39 - j0.90 


-0.03 + j0.12 


6960.0 


-0.05 H- j0.03 


-0.98 + j0.16 


-0.08 + j0.07 


72S0.0 


-0.19 + j0.03 


-0.08 + jO.97 


-0.18 + j0.15 


7600.0 


-0.13 + j0.14 


0.92 + jO.36 


-0.05 + j0.18 


7920.0 


0.01 + j0.07 


0.57 - j0.81 


0.01 + jO.ll 


8240.0 


-0.07 + jO.OO 


-0.63 - jO.76 


0.02 + jO.ll 


S560.0 


-0.12 + j0.04 


-0.90 + jO.42 


0.03 - jO.Ol 


88S0.0 


-0.10 + j0.07 


0.18 -t- jO.97 


-0.11 - j0.03 


9200.0 


0.00 + j0.19 


0.96 + j0.06 


-0.03 + j0.20 


9520.0 


0.10 + j0.14 


0.33 - j0.91 


0.11 + jO.ll 


9840.0 


0.11 - j0.04 


-0.81 - jO.55 


-0.02 - j0.09 


10160.0 


0.12 - j0.19 


-0.74 + jO.66 


O 

0 

1 

o 

b 


10480.0 


-0.03 - j0.13 


0.42 + jO.88 


0.03 + j0.02 


10800.0 


-0.02 - j0.14 


0.96 - j0.18 


-0.02 - j0.09 


11120.0 


-0.18 - j0.12 


0.06 - jO.98 


-0.16 + j0.03 


1 1440.0 


-0.14 - j0.03 


-0.91 - jO.33 


-0.08 + j0.02 


11760.0 


-0.02 - j0.14 


-0.57 + jO.79 


-0.04 - j0.04 


12080.0 


-0.04 - j0.04 


0.63 + jO.75 


0.04 + j0.06 


12400.0 


-0.02 - jO.lO 


0.86 - jO.45 


-0.02 - j0.04 
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Table 17. THE S - PARAMETERS OF FILTER + FIXTURE 



Freq.(MHz) 


5,. 




5:: 


6000.0 


-0.79 - jO.45 


0.00 4- jO.OO 


0.54 4 - jO.75 


6320.0 


0.89 - j0.14 


0.00 + jO.OO 


-0.36 - j0.S2 


6640.0 


-0.65 H- jO.62 


0.01 + jO.Ol 


-0.02 4 - j0.90 


6960.0 


0.11 - j0.91 


0.01 + jO.Ol 


0.05 - jO.S9 


72S0.0 


0.62 -h jo. 64 


0.01 -1- jO.Ol 


-0.24 4 - j0.90 


7600.0 


-0.80 + jO.45 


-0 02 4- j0.02 


0.02 - j0.91 


7920.0 


0.11 - jO.24 


-0.30 + jO.86 


0.00 4- jO.25 


8240.0 


-0.13 + j0.20 


-0.34 + jO.86 


-0.01 - j0.15 


8560.0 


-O.Oi jO.36 


-0.82 4- jO.24 


0.15 + jO.27 


8880.0 


0.0 1 j0.04 


-0.61 - jO.69 


-0.01 - j0.05 


9200.0 


0.31 + jO.OI 


0.56 - jO.65 


0.03 + j0.31 


9520.0 


0.27 jO.29 


0.76 4- j0.41 


-0.28 - j0.06 


9840.0 


-0.16 + jO.22 


-0.48 + jO.76 


0.07 - jO.23 


10160.0 


-0.06 + jO.23 


-0.73 - j0.50 


0.00 4- jO.29 


10480.0 


-0.06 - jO.26 


0.47 - jO.76 


-0.26 - j0.03 


10800.0 


0.21 - j0.02 


0.79 + jO.46 


-0.05 - j0.19 


11120.0 


-0.17 - j0.07 


-0.45 4- jO.78 


-0.14 - j0.02 


1 1440.0 


-0.14 + jO.OO 


-0.83 - jO.37 


0.14 + jO.Ol 


11760.0 


0.19 + j0.07 


0.16 - jO.86 


0.17 4- j0.05 


12080.0 


0.25 - jO.24 


0.72 - jO.28 


-0.48 - j0.13 


12400.0 


-0.15 + j0.70 


0.74 4- j0.18 


0.59 4- j0.12 
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Table 18 . THE S - PARAMETERS OF FILTER WITHOUT FIX- 
TURES IN DB UNITS 



Freq.(MHz) 








5:2 


6000.0 


1.57 


83.91 


88.79 


0.62 


6320.0 


1.96 


91.38 


105.37 


1.69 


6640.0 


0.47 


75.81 


78.12 


0.19 


6960.0 


0.87 


69.94 


69.39 


0.85 


72S0.0 


0.77 


54.40 


51.94 


0.86 


7600.0 


0.70 


29.99 


31.60 


1.12 


7920.0 


10.17 


0.57 


0.50 


9.15 


S 240.0 


11.32 


0.15 


0.16 


12.00 


8560.0 


8.62 


0.45 


0.45 


10.05 


SSSO.O 


15.77 


0.20 


0.08 


13.81 


9200.0 


9.66 


0.40 


0.20 


9.88 


9520.0 


10.97 


0.25 


0.33 


11.81 


9S40.0 


9.65 


0.27 


0.66 


9.11 


10160.0 


8.39 


0.45 


0.40 


10.81 


10480.0 


12.78 


0.15 


0.20 


12.01 


10800.0 


11.35 


0.43 


0.45 


11.87 


11120.0 


13.64 


0.30 


0.20 


14.39 


11440.0 


10.28 


0.31 


0.46 


16.09 


11760.0 


15.44 


0.60 


0.41 


15.84 


12080.0 


7.35 


1.40 


1.10 


9.98 


12400.0 


8.85 


1.45 


0.96 


11.51 
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Table 19. THE DE-EMBEDDED OUTPUT OF FILTER UNDER 
TEST 



Freq.(MHz) 


•^11 






5:, 


6000.0 


1.58 


37.03 


90.64 


0.28 


6320.0 


2.03 


41.18 


91.12 


1.15 


6640.0 


1.11 


41.62 


76.82 


0.14 


6960.0 


2.32 


59.94 


71.12 


0.90 


72S0.0 


0.68 


51.99 


53.34 


1.21 


7600.0 


0.21 


30.06 


31.45 


0.92 


7920.0 


10.50 


1.27 


1.51 


9.64 


8240.0 


10.99 


0.59 


0.60 


13.25 


8560.0 


7.97 


1.25 


1.18 


10.17 


8880.0 


12.83 


1.12 


0.55 


10.50 


9200.0 


9.57 


1.36 


0.70 


10.77 


9520.0 


8.80 


0.35 


0.94 


9.70 


9840.0 


10.83 


0.11 


0.81 


8.38 


10160.0 


7.58 


0.74 


0.90 


10.05 


10480.0 


11.96 


0.62 


0.75 


9.76 


10800.0 


12.39 


0.67 


0.73 


9.99 


11120.0 


11.49 


0.14 


0.20 


9.78 


1 1440.0 


8.94 


0.38 


0.76 


15.48 


11760.0 


13.73 


0.76 


0.94 


14.58 


12080.0 


8.20 


2.11 


1.41 


6.85 


12400.0 


8.41 


1.77 


1.38 


7.79 
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APPENDIX E. FORTRAN PROGRAM FOR DE-EMBEDDING THE 51 

OHM RESISTOR 



DIMENSION X(200),MDB11(200) ,PSE11(200) 
MDB12(200),PSE12(200) 
iMDB21(200),PSE21(200) 

MDB22(200) ,PSE22(200) ,MDB222(200) 

REAL MDB 1 1 , MDB 1 2 , MDB 2 1 , MDB2 2 , MDB 2 2 2 
MAGI 1 , MAG 12 , MAG2 1 , MAG22 
IM11,IM12,IM21,IM22 
COMPLEX S11,S12,S21,S22,A12A21,B12B21 

A11,A12,A21,A22,AA121,A12B12,A21B21 
B11,B12,B21,B22,BB121 
M11,M12,M21,M22 
E,FF,G,H,N,AABB 
PI=4.*ATAN(1. ) 

DO 10 1=1,21 

F=18. 5*10. **8+1. 5*10. **8*FL0AT( I) 
X(I)=F/1000000. 

CALL FIX1(A11,A12,A21,A22,A12A21,AMAG12,PSEA12,F) 
CALL TMFS(M11,M12,M21,M22,F) 

CALL FIX2(B11,B12,B21,B22,B12B21,BMAG12,PSEB12,F) 



TMAG=AMAG12*BMAG12 

TPSE=PSEA12+PSEB12 

XRE=TMAG*COS(TPSE) 

XIM=TMAG*SIN(TPSE) 

A12B12=CMPLX(XRE,XIM) 

A21B21=A12B12 

E=M12/(A12B12) 

FF=M21/(A21B21) 

G=(M11-A11)/A12A21 

H=(M22-B22)/B12B21 

N=( 1. +A22*G)*( 1. +B11*H)-A22*B11*E*FF 



S 1 1=( G*( 1 . +B 1 1*H ) -B 1 1*E*FF ) / N 

S12=E/N 

S21=FF/N 

S22=(H*(1.+A22*G)-A22*E*FF)/N 



QIM1=AIMAG(S11) 



QRE1=REAL(S11) 

QIM2=AIMAG(SI2) 

QRE2=REAL(S12) 

QIM3=AIMAG(S21) 

QRE3=REAL(S21) 

QIM4=AIMAG(S22) 

QRE4=REAL(S22) 

FIND A MAGNITUDE AND PHASE( ANGLE) 

MAG11=CABS(S11) 

MAG12=CABS(S12) 

MAG21=CABS(S21) 

MAG22=CABS(S22) 

MDB11(I)=ABS(20*ALOG10(MAG11)) 

MDB12(I)=ABS(20*ALOG10(MAG12)) 

MDB21(I)=ABS(20*ALOG10(MAG21)) 

MDB 2 2 ( I ) =AB S ( 2 0* ALOG 1 0 ( MAG2 2 ) ) 
PSE11(I)=ATAN2(QIM1,QRE1)*180. /PI 
PSE12(I)=ATAN2(QIM2,QRE2)*180. /PI 
PSE21(I)=ATAN2(QIM3,QRE3)*180. /PI 
PSE22(I)=ATAN2(QIM4,QRE4)*180. /PI 

WRITE(6,101) X(I),MDB11(I),MDB12(I),MDB21(I),MDB22(I), 
1 PSEll(I) ,PSE12(I),PSE21(I) ,PSE22(I) 

101 FORMAT(F6. 0,8( 1X,F7. 2)) 

10 CONTINUE 
END 



Q ■>W?'>W<’'3Wf'3VTV'>V’>V'>V‘>V’>V^Wf*'>V‘)V')V'>V^V'>V'>V‘>Wc^V‘>V'jWcVf‘5V*"3V’)V’A’ 

SUBROUTINE TMFS ( Dll ,D12 ,D21 ,D22 ,F) 

Q *Vf*>Wr>Wr>V7V’5Wf7V'5Wr'>V'>V^V'>V*>V'sV*>V*'>V'>Wr'>Wc'>V'>V'>Wc">Wc*3V*)V'>V 

COMPLEX D11,D12,D21,D22 
PI=4.*ATAN(1. ) 

INPUT DATA FILE IS DISK 9 
( BOTH FIXTURE + DEVICE ) 

READ(9,21)DB11,ANG11,DB21,ANG21,DB12,ANG12,DB22,ANG22 
21 FORMAT(8(F7. 2)) 

Vll=10.**(-DBll/20. ) 

V12=10.**(-DB12/20. ) 

V2 1=10. **(-DB2 1/20. ) 

V22=10.**(-dB22/20. ) 

RAD11=PI/180.*ANG11 
RAD12=PI/180. *ANG12 
RAD21=PI/180.*ANG21 
RAD22=PI/180.*ANG22 

QRE11=V11*C0S( RADII) 

QRE12=V12*COS(RAD12) 

QRE2 1=V2 l*COS( RAD2 1 ) 
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QRE22=V22*C0S(RAD22) 
QIM11=V11*SIN( RADII) 
QIM12=V12*SIN(RAD12) 
QIM21=V21*SIN(RAD21) 
QIM22=V22*SIN(RAD22) 

D11=CMPLX(QRE11,QIM11) 
D12=CMPLX(QRE12,QIM12) 
D2 1=CMPLX( QRE2 1 , QIM2 1 ) 
D22=CMPLX( QRE22 , QIM22 ) 

RETURN 

END 



C 

C 



Vc'>V'>V’5VVc’»V')V*')VV<’Vc'>V')V'>V’5V'A‘’5V')V'A’*')VVc’')V*Vc')V'5V*VcVr*'>V**’5V*')V')V'»V'5Vyf**Vc')VVc'>V'>V'jV*'>V'A’ 



SUBROUTINE FIX1(S11,S12 ,S21,S22,SS,SMAG12,SPSE12,F) 



*VcVfVfV-VrVc')V*V-'>VVf'>VV<’'5V')VVc')VVc'>VV-Vc’'jVV?')VVr’)VVr 



"j'oV'jWcVf'jWf’jWr'sWc’Vf'iVVf'jWf^VfVfjWc’Vc'A’ 



FOR FIXTURE A 

COMPLEX A(3,3),H(3,3) ,B(1,3),C(3),D(3,3),SS121 
COMPLEX S11,S12,S21,S22,SS 

CALL DATA1(A,B,F) 

CALL CPX(A) 

DO 10 1=1,3 

C(I)=A(1,I)*B(1,1)+A(2,I)*B(1,2)+A(3,I)*B(1,3) 
10 CONTINUE 

S11=C(1) 

S22=C(2) 

SS121=C(3) 

SS=S11*S22-SS121 

SIM=AIMAG(SS) 

SRE=REAL(SS) 

SMAGSS=CABS(SS) 

PSES S=ATAN2( S IM , SRE ) 

SPSE12=0. 5*PSESS 
SPSE21=SPSE12 
SMAG12=SMAGSS**0. 5 
SMAG21=SMAG12 

S12=(S11*S22-SS121)**0. 5 

S21=S12 

RETURN 

END 



C 

C 



Vr Vc Vc V? W Vc 



SUBROUTINE DATAlp,B,F) 

COMPLEX D(3,3) ,B( 1,3) ,RIN1,RIN2,RIN3,RL1,RL2,RL3 
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CALL RL0AD(RL1,RL2,RL3,F) 

CALL RIN(RIN1,RIN2,RIN3) 

D(l,l)=l. /RLl 

D(l,2)=l. /RL2 

D( 1,3) = 1. /RL3 

D(2,1)=RIN1 

D(2,2)=RIN2 

D(2,3)=RIN3 

D(3,l)=(-1. ,.0) 

D(3,2)=(-l. ,.0) 

D(3,3)=(-l. 0) 

B(1,1)=RIN1/RL1 

B(1,2)=RIN2/RL2 

B(1,3)=RIN3/RL3 

RETURN 

END 



Vr Vc * ")V *5 Wr V«"} V "} V "5 V ^ Wc”} Wrj W(T *5 V ^ V "jV -5 V 1 Wf * -jWf * *} Wf * 
SUBROUTINE RIN(RIN1 ,RIN2 ,RIN3) 

*/c^V'}WfVc'5Wc’*5W?Vf*V?’}V*)V'3V"5Wr'5Wr^Wr;V')WrVr5Wc-sWr5V 



COMPLEX RIN1,RIN2,RIN3 

CALL F1DAT(DB111,ANG111,DB112,ANG112,DB113,ANG113) 

PI=4. *ATAN( 1. ) 

INPUT DATA FILE IS DISKS 
(FIXTURE 1 ONLY) 

Vlll=10.**(-DBlll/20. ) 

V112=10. **(-dbi12/20. ) 

V113=10.**(-DB113/20. ) 

RADI 11=PI/ 180. *ANG111 
RADI 12=PI/ 180. *ANG112 
RADI 13=PI/ 180. *ANG113 

QRE111=V111*C0S(RAD111) 

QRE112=V112*C0S(RAD112) 

QRE113=V113*COS(RAD113) 

QIM111=V111*SIN(RAD111) 

QIM112=V112*SIN(RAD112) 

QIM113=V113*SIN(RAD113) 

RIN1=CMPLX(QRE111,QIM111) 

RIN2=CMPLX(QRE112,QIM112) 

RIN3=CMPLX(QRE113,QIM113) 

RETURN 

END 






■ “.V -:V ** Vr * Vf Vf ** -iWc *** Vf * * Vf **** 



SUBROUTINE F1DAT(DB111 , ANGl 1 1 , DB 1 12 , ANG112 ,DB113 , ANG113) 

y .y.y. ^ V' Vr Vr Vr Vc Vr V* Vr 'i: * V Vr Vr Vr Vr Vr ‘V Vr Vc Vc* V' Vc Vf Vf Vr V ' V* ixic V« Vr ^ 



INPUT DATA FILE IS DISK 8 
(For FIXTURE A) 

READ(8,20)DB111,ANG111,DB112,ANG112,DB113,ANG113 
FORMAT(6(F7. 2)) 

RETURN 

END 



y . y. ^ ^ ^ .t. y > y > y > y . y • y > y. y « y > Yc Vr Vc V* Vc 'V Vc Vv Vr Vr V' Vc Vr Vr Vr Vr ‘ V Vr Vc Vc Vr Vc 

SUBROUTINE FIX2(S1 1 , S12 , S21 , S22 ,SSS ,SHAG12, SPSE12 ,F) 






t. y.y. y. y. y.y^y^ y^ y 



. 



y«y« 



FOR FIXTURE 2 

COMPLEX A(3,3),H(3,3) ,B(1,3),C(3) ,D(3,3) ,SS121 
COMPLEX S11,S12,S21,S22,SSS 

CALL DATA2(A,B,F) 

CALL CPX(A) 

DO 10 1=1,3 

C(I)=A(1,I)*B(1,1)+A(2,I)*B(1,2)+A(3,I)*B(1,3) 

CONTINUE 

S11=C(1) 

S22=C(2) 

SS121=C(3) 

SSS=S11*S22-SS121 

SIM=AIMAG(SSS) 

SRE=REAL(SSS) 

SMAGSS=CABS(SSS) 

PSESS=ATAN2(SIM,SRE) 

SPSE12=0. 5*PSESS 
SPSE21=SPSE12 
SMAG12=SMAGSS**0. 5 
SMAG21=SMAG12 

S12=(S11*S22-SS121)**0. 5 

S21=S12 

RETURN 

END 



****yf*****VfVf*Vf*‘**'sWfV<’Vf** 

SUBROUTINE DATA2(D,B,F) 

COMPLEX D(3,3) ,B( 1,3) ,RIN1,RIN2,RIN3,RL1,RL2,RL3 
CALL RL0AD(RL1,RL2,RL3,F) 

CALL R2IN(RIN1,RIN2,RIN3) 



D(l,l) = l. /RLl 
D(l,2) = l. /RL2 
D(l,3)=l. /RL3 
D(2,1)=RIN1 
D(2,2)=RIN2 
D(2,3)=RIN3 
D(3,l)=(-1. ,.0) 

D(3,2)=(-l. ,.0) 

D(3,3)=(-l. 0) 

B(1,1)=RIN1/RL1 

B( 1,2)=RIN2/RL2 

B(1,3)=RIN3/RL3 

RETURN 

END 



C 

C 



V* VfVv y* V? y* y? y% y» y» y?y» y* y» y» y» y» y» y? y? 

SUBROUTINE R2IN(RIN1,RIN2,RIN3) 
y? y» y? y» y » y» y» y* y» y? y» y? * v y% y » y» y r y * y * y * y? y* y^ y ? y? *v 



COMPLEX RIN1,RIN2,RIN3 

CALL F2DAT(DB111,ANG111,DB112,ANG112,DB113,ANG113) 

PI=4.*ATAN(1. ) 

VI 11=10. **(-DB 11 1/20. ) 

V112=10.**(-DB112/20. ) 

V113=10. **( -DB113/20. ) 

RAD111=PI/180.*ANG111 
RAD112=PI/180.*ANG112 
RAD113=PI/180. *ANG113 

QRE111=V111*C0S(RAD111) 

QRE112=V112*COS(RAD112) 

QRE113=V113*COS(RAD113) 

QIM111=V111*SIN(RAD111) 

QIM112=V112*SIN(RAD112) 

QIM113=V113*SIN(RAD113) 

RIN1=CMPLX(QRE111,QIM111) 

RIN2=CMPLX(QRE112,QIM112) 

RIN3=CMPLX(QRE113,QIM113) 

RETURN 

END 



Q y? y? yfyfVfyr yry^y? yryfyryfyryryfyryf y?yfy»yf y? y? ycy? y? yryrvr yryrye-yf y? y^y? yryfyryryryf y^y? yryeyfyryryryey? yry? ycyp 

SUBROUTINE F2DAT(DB111 , ANGl 1 1 , DB 1 12 , ANG112 ,DB1 13 , ANGl 13) 



INPUT DATA FILE IS DISK 10 
(For FIXTURE B) 

READ( 10,20)DB111,ANG111,DB112,ANG112,DB113,ANG113 
20 FORMAT(6(F7. 2)) 
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RETURN 

END 



Vr*****'5V')V*'>VVcVfVr**'>ViV 

SUBROUTINE CPX(A) 

Vc “5 V * "5 V V? -5 V Vf ■5 V "5 V -sV "5 V ■5 V Vc * Vf it it 



INTEGER IP(3) 

COMPLEX A(3,3) ,W(3) ,DET(2) 

CALL CGEFA(A,3,3,IP,INFO) 

CALL CGEDI(A,3,3,IP,DET,W,01) 

RETURN 

END 



it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it 



SUBROUTINE RL0AD(RL1,RL2,RL3,F) 



COMPLEX RLl , RL2 , RL3 , ZO , ZLl , ZL2 , ZL3 
REAL L1,L2,L3 
PI=4. *ATAN(1. ) 

C=3. *10**8 
Z0S=68. 2 
ERE=2. 833 
ERES=(ERE)**0. 5 
L1=0. 04 
L2=0. 03 
L3=0. 02 

A1=C0S(2. *PI*F*L1*ERES/C) 

B1=SIN(2. *PI*F*L1*ERES/C) 

A2=C0S(2. *PI*F*L2*ERES/C) 

B2=SIN(2. *PI*F*L2*ERES/C) 

A3=C0S( 2. *PI*F*L3*ERES/C) 

B3=SIN(2. *PI*F*L3*ERES/C) 
R1=-Z0S*(A1/B1) 

R2=-Z0S*(A2/B2) 

R3=-ZOS*(A3/B3) 

ZL1=CMPLX(0. ,R1) 

ZL2=CMPLX(0. ,R2) 

ZL3=CMPLX(0. ,R3) 

Z0=CMPLX(50. ,0. ) 
RL1=(ZL1-Z0)/(ZL1+Z0) 
RL2=(ZL2-Z0)/(ZL2+Z0) 
RL3=(ZL3-Z0)/(ZL3+Z0) 

RETURN 

END 



Vc Vc Vr Vc V? Vr V? Vc Vf Vr Vc Vr Vc Vc Vr Vc Vc Vc Vc V? Vf 



THE INPUT DATA FILE FOR DISK 8 
3 DIFFERENT LOAD FOR FIXTURE A 



3. 07 


22. 7 


2. 28 


98. 2 


2. 14 


175. 1 


3. 82 


-63. 1 


2. 89 


29. 5 


2. 21 


114. 5 


3. 79 


-153. 9 


3. 77 


-49. 2 


2. 57 


55.4 


3. 46 


126. 4 


4. 12 


-138. 0 


3. 25 


-5. 1 


3. 20 


60. 5 


3. 75 


136. 6 


3. 86 


-74. 3 


3. 33 


-1. 7 


3. 31 


65. 1 


4. 44 


-156. 0 


4. 04 


-69. 9 


3. 40 


4. 9 


4. 42 


119. 8 


5. 39 


-152. 8 


3. 60 


-54. 2 


3. 98 


46. 6 


6. 17 


107. 6 


4. 55 


-119. 0 


3. 89 


-14. 2 


5. 32 


17. 1 


5. 78 


159. 8 


3. 89 


-66. 5 


4. 65 


-51. 0 


6. 39 


61. 8 


4. 15 


-120. 7 


3. 97 


-111. 9 


5. 76 


-25. 5 


5. 02 


175. 9 


4. 89 


-175. 9 


5. 70 


-89. 0 


8. 22 


99. 0 


7. 73 


104. 5 


4. 69 


-134. 2 


8. 95 


9. 5 


7. 36 


-17. 1 


3. 91 


173. 7 


6. 37 


-79.5 


6. 70 


-121. 4 


4. 48 


117. 4 


4. 35 


-147. 2 


4. 59 


170. 4 


6. 84 


47. 9 


3. 49 


157. 8 


4. 12 


114. 2 


8. 48 


-61. 4 


3. 37 


108. 4 


4. 88 


58. 9 


7. 57 


-174. 8 


3. 90 


58. 8 


7. 91 


-9. 1 


4. 92 


115. 1 


5. 55 


1. 4 


8. 42 


-137. 9 


3. 99 


62. 0 


6. 59 


-82. 6 



Vf Vf Vf V" Vf v» v» Vf V'f V* V* Vf Vf Vf V* ^ ^ Vf Vf Vf y* y«* y» y* y* y? y^ 

THE INPUT DATA FILE FOR DISK 9 
BOTH FIXTURE AND 51 OHM RES. 



11. 14 


-9. 3 


4. 60 


51. 0 


4. 40 


51. 4 


6. 45 


-150. 


11. 38 


-81. 2 


4. 71 


-1. 5 


4. 64 


0. 2 


6. 27 


173. 


10. 18 


-154. 3 


4. 80 


-53. 9 


4. 90 


-52. 2 


6. 05 


137. 


8. 47 


143. 1 


5. 35 


-106. 1 


5. 36 


-103. 3 


5. 86 


103. 


7. 13 


89. 7 


5. 84 


-152. 2 


5. 87 


-150. 5 


5. 51 


70. 


6. 35 


41. 0 


6. 25 


159. 0 


6. 30 


160. 2 


4. 99 


37. 


6. 21 


-5. 3 


6. 99 


112. 1 


6. 99 


113. 1 


4. 68 


3. 


6. 51 


-51. 8 


7. 49 


66. 6 


7.54 


68. 4 


4. 48 


-30. 


7. 33 


-102.0 


7. 25 


19. 4 


7. 32 


22. 0 


4. 43 


-65. 


8. 92 


-163. 3 


6. 65 


-26. 9 


7. 04 


-25. 9 


4. 94 


-101. 


11. 29 


112.5 


5. 81 


-79. 3 


7. 10 


-73. 4 


5. 75 


-133. 


12. 01 


12. 6 


6. 31 


-133. 3 


6. 46 


-121. 3 


6. 18 


-165. 


11. 59 


-61.5 


7. 25 


179. 3 


7. 30 


174. 6 


6. 91 


160. 


9. 15 


-93. 7 


6. 90 


135. 8 


6.45 


133. 3 


7. 83 


125. 


6. 27 


-142.5 


6. 56 


85. 4 


6. 35 


83. 2 


9. 07 


87. 


5. 53 


169. 5 


6. 25 


34. 5 


6. 00 


34. 1 


1. 78 


41. 


6. 07 


126. 3 


5. 85 


-17. 5 


5. 61 


-14. 6 


12. 72 


-22. 


7. 54 


87. 7 


5. 50 


-69. 5 


5. 30 


-64. 8 


12. 53 


-104. 


9. 91 


51. 8 


5. 39 


-121. 9 


5. 34 


-118. 0 


10. 00 


-174. 


10. 46 


23. 7 


5. 66 


-176. 7 


5. 70 


-171. 1 


7. 82 


131. 


11. 95 


44. 3 


6. 05 


131.4 


6. 10 


134. 9 


6. 53 


87. 



2 

1 

6 

2 

3 

0 

0 

7 

6 

3 

5 

5 

7 

8 

1 

0 

9 

0 

9 

2 

5 
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•jV-jV it itititit ititititititititit itit ititititit ititit ’jViWrsV it 

THE INPUT DATA FILE FOR DISK 10 
3 DIFFERENT LOAD FOR FIXTURE B 

it it it it it it it it itit it it it itit it it it it it it it it it it it it it it it it 



1. 97 


-117. 8 


2. 


32 


-26. 


9 


1. 41 


78. 


1. 89 


-176. 2 


2. 


46 


-91. 


9 


1. 71 


28. 


2. 00 


130. 2 


2. 


44 


-155. 


5 


2. 16 


-26. 


2. 27 


81. 9 


2. 


51 


147. 


3 


2. 48 


-84. 


2. 40 


34. 8 


2. 


41 


96. 


6 


2. 69 


-145. 


2. 62 


-13. 9 


2. 


25 


50. 


1 


2. 85 


152. 


3. 04 


-67. 3 


2. 


35 


6. 


9 


2. 85 


94. 


3. 87 


-130. 3 


2. 


28 


-35. 


9 


2. 75 


43. 


5. 05 


151. 6 


2. 


51 


-78. 


5 


2. 77 


-2. 


5. 41 


64. 0 


2. 


69 


-123. 


9 


2. 68 


-44. 


4. 91 


-11. 7 


3. 


25 


-175. 


1 


2. 68 


-82. 


3. 87 


-67. 6 


4. 


27 


126. 


2 


1. 93 


-121. 


2. 90 


-118. 9 


5. 


24 


56. 


1 


1. 91 


-162. 


2. 55 


-166. 3 


5. 


38 


-20. 


6 


2. 02 


157. 


2. 68 


147. 0 


4. 


49 


-90. 


2 


2. 42 


113. 


3. 33 


96. 3 


3. 


49 


-147. 


9 


3. 21 


64. 


4. 95 


36.4 


2. 


76 


163. 


0 


4. 49 


7. 


7. 33 


-43. 4 


2. 


43 


118. 


9 


5. 70 


-60. 


7, 51 


-139. 7 


2. 


36 


77. 


5 


5. 80 


-134. 


5. 35 


143. 6 


2. 


55 


36. 


2 


4. 89 


158. 


3. 67 


87. 9 


3. 


32 


-8. 


6 


3. 72 


105. 



5 

1 

2 

3 

9 

4 

8 

2 

9 

3 

4 

1 

1 

0 

2 

6 

5 

8 

9 

8 

1 



SO 



APPENDIX F. FORTRAN PROGRAM FOR DE-EMBEDDING THE 

FILTER 



DIMENSION X(21),MDB11(21),PSE11(21) 

MDB12(21) ,PSE12(21) 

NDB2K21) ,PSE21(21) 

MDB22(21) ,PSE22(21) 

REAL MDB11,NDB12,MDB21,MDB22 
MAG 1 1 , MAG 12 , MAG2 1 , MAG22 
IM11,IM12,IM21,IM22 
COMPLEX S11,S12,S21,S22 
A11,A12,A21,A22 
B11,B12,B21,B22,W 
D11,D12,D21,D22 
C11,C12,C21,C22 
F11,F12,F21,F22 
M11,M12,M21,M22 
H11,H12,H21,H22 

PI=4.*ATAN(1. ) 

DO 10 1=1,21 

INPUT DATA FILE IS DISKS 
(FIXTURE 1 ONLY) 

READ(8,20)DB11,DB12,DB21,DB22,ANG11,ANG12,ANG21,ANG22 
20 FORMAT(8(F7. 2)) 

Vll=10.**(-DBll/20, ) 

V12=10,**(-DB12/20. ) 

V21=10.**(-DB21/20. ) 

V22=10.**(-DB22/20. ) 

RAD11=PI/180.*ANG11 

RAD12=PI/180.*ANG12 

RAD21=PI/180.*ANG21 

RAD22=PI/180.*ANG22 

QRE 1 1= V 1 1*C0 S ( RAD 1 1 ) 

QRE12=V12*C0S(RAD12) 

QRE21=V21*C0S(RAD21) 

QRE22=V22*C0S(RAD22) 

QIM11=V11*SIN(RAD11) 

QIM12=V12*SIN(RAD12) 

QIM21=V21*SIN(RAD21) 

QIM22=V22*SIN(RAD22) 

A11=CMPLX(QRE11,QIM11) 

A12=CMPLX(QRE12,QIM12) 



81 



A2 1=CMPLX( QRE 2 1 , QIM2 1 ) 

A22=CMPLX(QR£22,QIM22) 

C "S"PARAMETERS TO "T"PARAMETERS 

B11=-(A11*A22-A12*A21)/A21 
B12= A11/A21 
B21=-A22/A21 
B22=l/A21 

C INVERSE OF FIXTURE 1 "T" PARAMETER 

W= 1/(B11*B22-B21*B12) 

C11=W*B22 
C12=-W*B21 
C21= -W*B12 
C22= W*B11 

CALL OF 'T' PARAMETERS OF BOTH FIXTURE + FILTER 

CALL TMFS(D11,D12,D21,D22,I) 

C (INVERSE OF FIXTURE 1) * (BOTH FIXTURE+DEVICE) 

H11=C11*D11+C21*D12 

H12=C12*D11+C22*D12 

H21=C11*D21+C21*D22 

H22=C12*D21+C22*D22 

CALL OF INVERSE 't' PARAMETER OF FIXTURE 2 
CALL FIX2(F11,F12,F21,F22) 

C CASCADE OF FIXTURE 1, FILTER AND FIXTURE 2 



M11=H11*F11+H21*F12 

M12=H12*F11+H22*F12 

M21=H11*F21+H21*F22 

M22=H12*F21+H22*F22 



C "T"PARAMETERS TO "S"PARAMETERS 
S11=M12/M22 

S12= M11-(M12*M21)/M22 
S21= 1/M22 
S22= -M21/M22 

QIM1=AIMAG(S11) 

QRE1=REAL(S11) 

QIM2=AIMAG(S12) 

QRE2=REAL(S12) 

QIM3=AIMAG(S21) 

QRE3=REAL(S21) 

QIM4=AIMAG(S22) 

QRE4=REAL(S22) 
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c 



DE -EMBEDDED OUTPUT 



MAG11=CABS(S11) 

MAG12=CABS(S12) 

MAG21=CABS(S21) 

MAG22=CABS(S22) 

MDB11(I)=ABS(20*AL0G10(MAG11)) 

MDB12( I)=ABS( 20*ALOG10(MAG12) ) 
MDE21(I)=ABS(20*AL0G10(MAG21)) 
MD322(I)=ABS(20*AL0G10(HAG22)) 
PSE11(I)=ATAN2(QIM1,QRE1)*180. /PI 
PSE12(I)=ATAN2(QIM2,QRE2)*180. /PI 
PSE21(I)=ATAN2(QIM3,QRE3)*180. /PI 
PSE22(I)=ATAN2(QIM4,QRE4)*180. /PI 

N=5680 

X(I)=N+320. *FLOAT(I) 

WRITE(6, 100) X(I) ,MDB11( I) ,MDB12(I) ,MDB21( I) ,MDB22( I) , 
1 PSEll(I) ,PSE12(I),PSE21(I) ,PSE22(I) 

100 FORMATC ' ' ,2X,F6. 0,8(1X,F7. 2)) 

10 CONTINUE 
END 



Q •jVVr'jVVr'jVVc'jVyr^-sV****************’************* 

SUBROUTINE TMFS (Dll ,D12 ,D21,D22 , I) 

Q •5V'>V';V')V':V'5VVf'>V'5Vyr'5V';V')V*')V')V*yf'>V'>V'>V‘)V';VyrV*-'>Vyc‘jV')V*V'C'>Vyr')Vyr'5Vypy«’yr 

COMPLEX A11,A12,A21,A22 
COMPLEX D11,D12,D21,D22 
PI=4, *ATAN( 1. ) 

INPUT DATA FILE IS DISK 9 
(BOTH FIXTURE + DEVICE ) 

READ(9,21)DB11,DB12,DB21,DB22,ANG11,ANG12,ANG21,ANG22 
21 FORMAT(8(F7. 2)) 



Vll=10. **(-DB11/20. ) 
V12=10.**(-DB12/20. ) 
V2 1=10. **(-DB2 1/20. ) 
V22=10.**(-DB22/20. ) 
RAD11=PI/180.*ANG11 
RAD12=PI/180.*ANG12 
RAD21=PI/180.*ANG21 
RAD22=PI/180. *ANG22 

QRE11=V11*C0S(RAD11) 
QRE12=V12*C0S(RAD12) 
QRE2 1=V2 1*C0S ( RAD2 1 ) 
QRE22=V22*C0S( RAD22 ) 
QIM11=V11*SIN(RAD11) 
QIM12=V12*SIN(RAD12) 
QIM21=V21*SIN(RAD21) 
QIM22=V22*SIN(RAD22) 
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A11=CMPLX(QRE11,QIM11) 

A12=CMPLX(QRE12,QIM12) 

A2 1=CMPLX( QRE21 , QIM2 1 ) 
A22=CMPLX( QRE22 , QIM22 ) 

D11=-(A11*A22-A12*A21)/A21 

D12= A11/A21 

D21=-A22/A21 

D22=l/A21 

RETURN 

END 



*^V';V*;V*‘>V**')V^V'}V‘5V‘>VV(r'>V*'iV‘)V’3V‘5V‘>V**>VVe‘)V')V')V‘5V‘5V‘jV‘5V* 

SUBROUTINE FIX2(F11 ,F12 ,F21 ,F22) 

Vf'>V';V*VcV-*^V*)Vy?*‘V‘V'jVVrVfyc’'jVy?’jV‘>V*5VVc’‘3Wf’)Wr‘jV‘;V'5V‘>Wc 



COMPLEX A11,A12,A21,A22 
COMPLEX B11,B12,B21,B22,W 
COMPLEX C11,C12,C21,C22 
COMPLEX F11,F12,F21,F22 
PI=4, *ATAN( 1. ) 

INPUT DATA FILE IS DISK 10 
(FIXTURE 2 ONLY) 

READ(10,20)DB11,DB12,DB21,DB22,ANG11,ANG12,ANG21,ANG22 
FORMAT(8(F7. 2)) 

Vll=10.**(-DBll/20. ) 

V12=10. **(-DB 12/20. ) 

V21=10.**(-DB21/20. ) 

V22=10. **( -DB22/20. ) 

RAD11=PI/180.*ANG11 
RAD12=PI/180. *aNG12 
RAD21=PI/180.*ANG21 
RAD22=PI/180.*ANG22 

QRE11=V11*C0S(RAD11) 

QRE 1 2=V 1 2*C OS ( RAD 1 2 ) 

QRE 2 1=V2 l*COS( RAD2 1 ) 

QRE22=V22*COS( RAD22 ) 

QIM11=V11*SIN(RAD11) 

QIM12=V12*SIN(RAD12) 

QIM21=V21*SIN(RAD21) 

QIM22=V22*SIN(RAD22) 

A11=CMPLX(QRE11,QIM11) 

A12=CMPLX(QRE12,QIM12) 

A2 1=CMPLX( QRE2 1 , QIM2 1 ) 

A22=CMPLX(QRE22,QIM22) 



B 1 1=- ( A1 1*A22 -A12*A2 1 ) / A2 1 
B12= A11/A21 



B21=-A22/A21 

B22=l/A21 

W= 1/(B11*B22-B21*B12) 

F11=W*B22 

F12=-W*B21 

F21= -W*B12 

F22= W*B11 

RETURN 

END 



Q VnWc”sV*'3Wc'>Wc**’>V'jV’3V'>Wc'5V*'3V’A’'3Wf ***■}'?**** 

c The Input Data File for Disk.8 

r '5V’sV‘}V*Vc’*.V*/V')Wf}V';W?')V'>V*'>V'>W-'>Wc**5Wc'5V'>V'>V')Wr*’}V 



14. 48 


. 19 


0. 


10 


16. 22 


116. 0 


34. 2 


34. 5 


125. 


19. 83 


. 20 


0. 


11 


14. 01 


154. 9 


-38. 8 


-39. 5 


148. 


22. 21 


. 15 


0. 


11 


18. 46 


146. 2 


-113. 3 


-114. 3 


105. 


24. 81 


. 10 


0. 


14 


19. 94 


142. 8 


170. 8 


170. 2 


139. 


14. 41 


. 20 


0. 


15 


12. 81 


171. 7 


94. 6 


95. 2 


140. 


14. 49 


. 09 


0. 


11 


14. 74 


133. 1 


21. 3 


21. 0 


104. 


22. 80 


. 05 


0. 


09 


19. 44 


80. 8 


-54. 8 


-54. 3 


82. 


23. 32 


. 10 


0. 


10 


19. 30 


179. 7 


-129. 7 


-129. 3 


78. 


17. 71 


. 10 


0. 


15 


30. 65 


160. 7 


154. 9 


155. 7 


-16. 


18. 60 


. 15 


0. 


11 


18. 74 


144. 4 


79. 2 


79. 8 


-162. 


14. 23 


. 30 


0. 


30 


13. 76 


88. 8 


3. 7 


4. 5 


99. 


15. 56 


. 26 


0. 


30 


16. 14 


54. 1 


-70. 3 


-69. 4 


43. 


19. 00 


. 17 


0. 


27 


20. 31 


-19. 0 


-145. 9 


-144. 3 


-102. 


12. 90 


. 10 
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